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The  Electromagnetic  Pulse  (EMP)  produced  by  a  high- 
altitude  nuclear  blast  presents  a  severe  threat  to  electronic 
systems  due  to  its  extreme  characteristics.  To  test  the 
vulnerability  of  large  systems,  such  as  airplanes,  missies,  or 
satellites,  they  must  be  subjected  to  a  simulated  EMP  en¬ 
vironment. 

One  type  of  simulator  that  has  been  used  to  approxi¬ 
mate  the  EMP  environment  is  the  Large  Parallel-Plate 
Bounded-Wave  Simulator.  It  is  a  guided  wave  simulator 
which  has  the  properties  of  a  transmission  line  and  supports 
a  single  TEM  mode  at  sufficiently  low  frequencies.  This 
type  of  simulator  consists  of  finite-width  parallel-plate 
waveguides,  which  are  excited  by  a  wave  launcher  and  ter¬ 
minated  by  a  wave  receptor. 

This  study  addresses  the  field  distribution  within  a  finite- 
width  parallel-plate  waveguide  that  is  matched  to  a  conical 
tapered  waveguide  at  either  end.  Characteristics  of  a  parallel- 
plate  bounded-wave  EMP  simulator  were  developed  using 
scattering  theory,  thin-wire  mesh  approximation  of  the 
conducting  surfaces,  and  the  Numerical  Electronics  Code 
(NEC).  Background  is  provided  for  readers  to  use  the  NEC 
as  a  tool  in  solving  thin-wire  scattering  problems. 
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CHAPTER  1 


INTRODUC  TION  R)  PARALLEL-PI. A  I  E  EMP  SIMULATORS 


I  he  Electromagnetic  Pulse  'EMP)  produced  by  a  high-altitude  nuclear  burst 


presents  a  severe  'hreat  to  electronic  systems  due  to  its  extreme  characteristics.  In 


order  to  test  the  vulnerability  of  large  systems,  such  as  airplanes,  missiles,  or 


satellites,  they  must  be  subjected  to  a  simulated  EMP  environment,  which  is 


characterized  by  a  high  peaked  transient  pulse  with  a  planar  wavefront  that  will 


uniformly  illuminate  a  test  object.  One  type  of  simulator  that  has  been  used  to 


approximate  the  EMP  environment  is  the  Large  Parallel-Plate  Bounded-Wave 


Simulator.  The  parallel-plate  bounded-wave  simulator  is  a  guided  wave  simulator 


which  has  the  properties  of  a  transmission  line  and  supports  a  single  TEM  mode  at 


sufficiently  low  frequencies.  Ideally,  it  has  a  uniform  characteristic  impedance  that  is 


independent  of  frequency  along  the  direction  of  propagation,  and  can  be  terminated  by 


a  matched  impedance  providing  a  pure  traveling  wave. 


Parallel-plate  bounded-wave  EMP  simulators  consist  of  finite-width  parallel- 


plate  waveguides,  which  are  excited  by  a  wave  launcher  and  terminated  by  a  wave 


receptor.  Typically,  conical  tapered  waveguides  are  matched  to  the  parallel-plate 


region  in  order  to  excite  and  terminate  a  plane  wave  in  the  parallel-plate  region.  This 


geometry  is  illustrated  in  Figure  1.1  (all  tables  and  figures  are  located  at  the  end  ol 


each  chapter).  The  apex  of  the  wave  launcher  is  excited  by  a  transient  source  which 


has  a  very  fast  rise  time  that  is  significantly  smaller  than  the  transit  time  of  the 


waveguide,  on  the  order  of  1-5  nanseconds.  and  a  fall  time  on  the  order  of  I 


microsecond.  In  theoretical  computations,  the  pulse  is  usually  represented  as  a  double 


exponential  waveform  and  has  broadband  characteristics. 


The  problem  being  considered  in  this  thesis  is  the  solution  of  the  field  distribution 


within  a  finite-width  parallel-plate  waveguide  that  is  matched  to  a  conical  tapered 


wasias 


waveguide  at  either  end.  The  complete  field  representation  is  a  superposition  of 
discrete  and  continuous  modes.  At  very  low  |  requencies.  where  the  plate  separation  is 
much  less  than  one-half  of  a  wavelength,  only  the  I  EM  modes  ot  the  discrete 
spectrum  will  propagate  along  the  waveguide.  By  assuming  a  finite-width  parallel- 
plate  waveguide  of  infinite  length,  the  TEM  field  distribution  can  be  obtained  by 
means  of  a  quasi-static  approximation  (1.2.3).  This  analogy  is  usel  ul  in  establishing 
the  characteristic  impedance  and  the  field  distribution  of  the  line  for  low  frequencies. 
Furthermore,  quasi-static  analysis  on  the  conical  tapered  waveguide  has  been 
performed  (4.5)  by  use  of  stereographic-pro jection  and  conformal-mapping  techniques. 

Under  a  transient  pulse  excitation  of  the  parallel-plate  simulator,  there  will  be 
significant  energy  distributed  throughout  a  very  broad  band,  extending  into  the  high 
frequency  region.  Unfortunately,  for  all  practical  purposes,  the  plate  separation  is 
quite  large  and  the  bandwidth  of  single-mode  operation  is  much  smaller  than  that  of 
interest  under  pulse  exalaiton.  As  a  result,  higher-order  TE  and  TM  modes  will 
propagate  along  the  waveguide.  Furthermore,  in  the  high-frequency  region,  the  open 
simulator  will  act  as  a  radiator  and  will  no  longer  support  purely  guided  modes. 

In  order  to  develop  an  accurate  approximation  of  the  field  distribution  within  the 
simulator.  non-TEM  modes  must,  be  studied.  Initially,  assume  the  parallel-plate 
section  to  be  inhnitely  wide  and  long.  In  this  case,  the  modes  can  be  thought  of  as 
those  existing  in  a  rectangular  waveguide  whose  width  goes  to  infinity  (6).  For  lower 
:  requencies,  it  has  been  shown  that  a  detailed  field  mapping,  carried  out  at  a  single 
CW  !  requency.  can  be  well  approximated  with  the  superposition  of  propagating  modes 
when  using  the  rectangular  waveguide  approximation.  However,  if  further  accuracy 
is  desired,  a  more  complete  analysis  can  be  derived  by  computing  the  discrete  modes  o! 
,i  very  long,  tinite-width  parallel-plate  waveguide.  This  problem  was  first  approached 
by  assuming  either  narrow  plates  (7)  or  wide  plates  (8.9,10).  In  (7).  Mann  found  that 


for  narrow  plates  (width  <  <  separation),  the  TE  modes  suffer  greater  attenuation  as 
they  propagate  than  the  I'M  modes.  Furthermore,  it  was  found  that  for  wide  plates 
(width  >>  separation),  the  1M  modes  have  a  greater  attenuation  than  the  IE  modes 
(8  d.it)!.  Therefore,  for  the  ease  ol  an  arbitrary  width,  one  should  expect  the  existence 
of  both  propagating  TE  and  TM  modes  inside  the  parallel-plate  region. 

A  more  complete  analysis  was  performed  at  the  University  of  Illinois  [l  1,12.  IJ] 
on  the  study  of  the  source  excitation  of  the  finite-width  parallel-plate  waveguide. 
This  analysis  involved  the  computation  of  the  total  E  field,  as  a  superposition  of  the 
discrete  and  continuous  spectrum.  The  results  of  [11.12]  are  valid  for  waveguides 
w'hose  cross  sections  are  large  compared  to  a  wavelength,  whereas  the  results  of  [lil 
are  valid  for  a  small-to-moderate  cross  section. 

It  is  important  to  understand  that  these  techniques  approximate  the  fields  for  an 
infinitely  long  parallel-plate  waveguide.  However,  the  parallel-plate  region  of  ’he 
simulator  is  finite  in  length  and  is  excited  and  terminated  by  conical  transmission 
lines.  Unfortunately,  the  previous  quasi-static  analyses  do  not  account  for  the 
generation  of  spurious  modes  due  to  the  mismatch  at  the  interface  of  the  two 
waveguides.  A  paper  written  by  John  Lam  [14]  develops  an  analytic  expression  for 
the  TEVI  modes  of  a  semi-infinite,  finite-width,  parallel-plate  waveguide,  excited  by  a 
conical  transmission  line.  Although  the  formulation  does  not  include  the  effects  of 
anv  reflected  fields,  it  does  yield  insight  into  the  excitation  process  of  the  parallel- 
plate  section.  It  should  be  realized  that  interfacing  the  two  waveguides  requires  the 
matching  of  a  spherical  wavefront,  which  has  £9  and  H $  components,  into  a  planar 
wavelront  supported  by  the  parallel-plate  region.  Along  the  vertical  plane  y  «  a 
constant,  the  I  idds  are  matched  across  the  transition  as  £.  =  —  E 9  sin  0,  £v  =  £9  cos  0. 
H t  =  -H $  cos  6  and  Hy  =  —  H $  sin  fr.  <b  represents  the  angle  between  the  center  axis 
and  the  observation  point.  0  its  elevation,  and  the  origin  is  at  the  apex  of  the  conical 


waveguide.  As  a  result,  there  will  exist  both  transverse  and  longitudinal  electric  and 
magnetic  fields  at  the  interface,  contributing  to  a  superposition  of  TEM.  TE.  and  TM 
modes.  However,  tor  gradual  tapered  end  sections  the  mismatch  is  relatively  small 
and  the  higher-order  modes  are  quite  small  for  low  frequencies. 

Furthermore,  at  the  transition,  there  exists  a  certain  time  dispersion  in  matching 
the  spherical  wave  into  a  plane  wave  [15]  since  the  cross  section  is  not  a  plane  of 
constant  phase.  Referring  to  Figure  1.2.  a  dispersion  distance  is  defined  and  is  related 
to  the  time  dispersion  as 

A  =  Ap  +  Ay  .  (1.1) 

'X  hen  the  wavelength  is  less  than  or  equal  to  four  times  the  dispersion  distance,  the 
higher-order  modes  produced  become  significant  in  amplitude  leading  to  a  large 
distortion  of  the  waveform.  Therefore,  the  dispersion  distance  determines  the  so- 
called  “band  width"  of  the  simulator. 

I  he  analyses  that  have  been  reviewed  thus  far  in  this  chapter  all  offer  excellent 
physical  and  theoretical  interpretations  of  the  parallel-plate  bounded-wave  EMP 
simulator  by  considering  ideal  situations.  Through  such  modal  analyses,  the 
explanation  of  notches  in  the  magnetic  field  has  been  determined  [16],  optimum  loads 
designed,  and  other  standards  in  design  have  been  developed  [17].  Furthermore, 
accurate  field  mappings  have  been  generated,  although  they  require  the  use  of 
experimental  data  for  their  generation.  This  is  quite  efficient  for  simulators  that  have 
been  constructed  and  for  frequencies  well  within  the  simulator  band-width.  When 
designing  a  parallel-plate  EMP  simulator,  the  engineer  has  much  more  to  consider  than 
i he  transmission  line  properties  of  the  simulator  itself.  It  is  very  important  to 
understand  that  once  a  test  object  is  placed  within  the  parallel-plate  region  the  relative 
amplitudes  of  'he  modes  may  differ  greatly  from  those  calculated  in  an  empty 


working  volume.  Therefore,  one  cannot  treat  the  calculated  fields  of  the  empty 
working  volume  as  the  incident  held  that  establishes  the  charge  and  current  on  the 
test  object.  This  so-called  "simulator-obstacle  interaction”  can  be  analysed  only  lor 
certain  geometries  [  1 8. 1  T20.2 1  j.  such  as  an  infinite  cylinder  in  a  parallel-plate 
waveguide.  Furthermore,  the  design  engineer  should  be  interested  in  the  tields  that 
wiil  be  radiated  into  the  external  environment  over  the  entire  frequency  spectrum. 
There  will  be  a  number  of  resonant  frequencies  that  exist  due  to  various  mismatches 
in  geometries  and  loads,  as  well  as  to  resonances  of  the  parallel-plate  region. 

For  the  design  of  an  alternate  geometry  or  configuration,  an  m-depth  analysis 
may  be  necessary  and  would  prove  to  be  quite  time-consuming  and  difficult. 
However,  to  develop  further  theoretical  insight  into  the  problem,  the  ability  to  easily 
perform  an  accurate  analysis  using  a  dynamic  field  solution  could  be  quite  useful. 
1  he  purpose  of  this  thesis  is  to  evaluate  the  well-known  Numerical  Flectromagnetics 
Code  (NEC),  developed  by  Lawrence  Livermore  Laboratory  (22).  as  a  potential  tool  in 
generating  a  readily  available  analysis  of  the  parallel-plate  bounded-wave  simulator 
in  the  frequency  domain.  The  NEC  code  is  a  user-oriented  code  that  is  widely 
distributed,  and  it  the  problem  is  interpreted  correctly,  it  can  be  a  very  useful  tool  for 
many  electromagnetic  scattering  and  antenna  problems. 

In  Chapter  2.  NEC  is  examined  and  the  validity  of  solving  scattering  problems 
using  'he  thin-wire  approximation  is  developed.  More  specifically,  the  two-wire 
mom  me  EMP  simulator,  as  well  as  a  t wo-dimensional  plate  scattering  problem,  are 
analx/ed.  The  emphasis  of  the  analogies  is  placed  on  interpreting  the  results  and 
appivmg  ’hem  to  the  development  of  a  wire  mesh  model  of  the  bounded-wave 


parallel-plate  I  MP  simulator. 

In  Chapter  L  'he  wire  mesh  approximation  ot  the  surface  of  the  conducting  plates 
.»t  the  parallel  plate  simulator  is  devised.  The  tieid  distribution  within  the  para  1  lei  - 


plate  region  due  to  a  continuous-wave  excitation  across  the  source  gap  is  then 
computed  using  NEC.  These  results  are  interpreted  and  the  limitations  of  the  wire 
mesh  approximation  are  investigated.  Finally,  various  mechanisms  contributing  to 
phenomena  such  as  the  simulator  band  width  and  simulator-obstacle  interaction  are 


CHAPTER  2 


SOLUTION  OF  THIN-WIRE  SC  A  UTERING  PROBLEMS  USING  THE 
NUMERICAL  ELECTROMAGNETICS  CODE 


This  chapter  is  intended  to  give  the  reader  a  useful  backround  to  employ  the 
Numerical  Electromagnetics  Code  (NEC)  as  a  too!  in  solving  thin-wire  scattering 
problems.  Initially,  a  brie!  theoretical  back  round  of  NEC  will  be  presented.  Then,  the 
validity  ot  its  solutions  will  be  analyzed  for  the  following  problems:  Thin  monopole 
over  a  PEC  ground  plane,  the  rhombic  EMP  simulator,  and  the  scattering  of  a  thin 
square  PEC  plate.  The  goal  of  this  analysis  is  to  determine  how  to  best  model  a 
conducting  surface  using  a  thin-wire  approximation,  as  well  as  how  to  best  model 
source  excitations.  With  the  understanding  of  how  to  interpret  the  results  of  the 
thin-wire  approximation,  this  analysis  will  provide  useful  insight  into  the  dynamic 
solution  of  the  parallel-plate  bounded-wave  EMP  simulator  using  NEC.  which  is 
presented  in  Chapter  i. 

2.1  I  heoretical  Background  of  NEC 

In  this  section.  NEC  will  be  used  to  solve  thin-wire  scattering  problems,  which 
involves  solving  for  the  induced  currents  on  the  wire's  surface  due  to  some  incident 
held.  However,  the  method  of  solution  used  by  NEC  must  be  understood  before  the 
•  tun -wire  scattering  problem  can  properly  be  developed.  The  complete  analysis  ot 
\[.(  , s  provided  in  (22).  The  following  presents  an  overview  of  the  theory  used  and 
on",  .lies  '  he  basis  I  or  our  analysis. 

The  Numerical  Electromagnetics  Code  is  based  on  a  moment  method  solution  ol 
•ne  Electric  Field  In'egral  Equation  (EFIE)  for  thin-wire  scatterers  and  the  Magnetic 
Field  Integral  Equation  (MF1E)  for  closed  volumes.  In  this  analyis.  only  the  solution 


of  the  EFIE  will  be  considered.  The  EFIE  used  by  NEC  is  based  on  the  thin-wi re 
approximation  in  which  the  surface  current  can  be*  assumed  to  be  translationailv 
invariant  and  a  1  unction  ol  the  direction  ot  its  wire  axis  only.  I  h is  is  valid  il  the 
radius  .a  'he  'Aire  scourer  is  much  smaller  than  a  wavelength.  I  hen  the  surluce 
current  can  he  approximated  by  a  single  ti lament  on  the  ware  axis,  and  the  boundary 
conditions  can  be  applied  on  the  ware  surface.  This  is  better  known  as  the  Extended 
Boundary  Condition  [2  <1.  Assuming  a  PEC'  surface,  the  applied  boundary  condition  on 
•he  tangential  1  ield  is  limited  to  the  axially  directed  field  and  is  expressed  as 


f:  +  r,s  =  o 


(2.1.1  .' 


where  F  is  the  axially  directed  incident  field  and  F 2  is  the  axially  directed  scattered 
1  ie Id  on  the  conducting  surface  which  is  produced  by  the  current  I  ilament  on  the  ware 
axis.  With  the  use  ol  ’he  Extended  Boundary  Condition,  the  EFIE  can  then  be  reduced 
to  a  lint-  integral,  integrating  over  the  ware  axis,  and  can  be  expressed  as 


- f/(z  )  +  k2  G 

l  dr- 


:  G  (  2  .2  )dz  +  /  (  2  tZ  (  2 


(2.1.2! 


w  here 


/  (  :  )  =  the  unknown  axial  current 


G(: 


,  - ;*  V a  1  *  1  :  a 
(’ 

4  jVj  :  +  ( :  —  j  )’ 


is  the  thin-wire  approximation  of  the  tree- 


space  t/reen's  !  unction  along  the  wire  axis 

/( i )  -  the  dist  ri billed  impedance  loading  along  the  wire  surt  ace 

a  =  the  wire  radius  and 
*>  — 

k  =  -1_  =  the  I  rtr-space  wave  number, 
x 

1  he  left-hand  side  of  Equation  <2.1.2)  is  the  axially  directed  incident  f  ie  Id  on  the 
ware  surtace  and  is  assumed  to  be  known  while  the  axial  current  is  the  unknown. 
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However,  this  cannot  be  solved  analytically,  in  general,  and  must  be  solved 
numerically.  The  numerical  solution  oi  the  EFIE  is  then  obtained  with  the  use  ol  the 
moment  method  solu'ion  (24). 

Inc  unknown  axial  current  is  expanded  into  a  truncated  series  ol  basis  I  unctions 
weighted  b\  unknown  constant  coet  I  icients.  The  basis  function  was  chosen  to  be  a 
weighted  three-term  basis  I  unction,  with  the  n:h  basis  function  represented  as 

/,(*)  =  .\,  +  U,  mildl+r,  las(d).  (2.1.3) 

There! ore.  three  unknowns  have  been  introduced  and  three  conditions  are  needed  to 
solve  tor  the  constant  weighting  coel  ficients. 

Initially .  the  w  ire  geometry  is  broken  up  into  a  set  of  N  individual  wire  segments, 
each  of  length  A,  .  Each  basis  function  has  a  finite  support  lying  on  an  interval  that  is 
centered  on  a  wire  segment,  and  extends  to  the  center  of  adjacent  segments  at  either 
end  ol  the  specified  wire  segment.  The  inner  product  of  both  sides  of  the  EFIE  is 
taken  with  a  set  ol  N  testing  functions.  The  testing  functions  are  chosen  to  be 
mathematical  point  functions,  where  then'*  point  function  is  located  at  the  center  of 
the  o  ’1  wire  segment.  The  integral  over  the  point  function  can  be  evaluated 
anaiv '  teal  I  y .  and  the  current  basis  functions  can  be  rewritten  as 

/,  (  r  )  =  .4,  +5,  sin  *  ( z  —2n  )  +  Cn  cos  k  ( z  — z,  )  (2.1.4) 

u  m  u  :  ;  -z,  i  <  -1  and  .%  is  the  center  of  the  n !h  segment.  Now.  two  of  the  three 

unknowns  can  be  satisfied  by  enforcing  charge  and  current  continuities  at  the 

A, 

endpoints  oi  •  ne  w,re  segments  (r  =  ±  —  ).  where  the  line  charge  density  is  related 


’o  \he  current  by  the  continuity  equation 


-  /  Uty 
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Iheretore.  with  the  use  of  the  three-term  basis  t  unc t ;on .  problems  such  as  open -ended 
and'or  'apered-wire  structures,  as  well  as  multiple  wire  nmc'.ons.  can  be  handled, 
wi'hin  limitations  oi  ' he  approximat ions. 

At  'he  junction  ot  two  or  more  wire  segments,  current  continurv  ,s  enforced  nv 
Kirchotf's  law.  which  states  that  the  sum  of  all  the  currents  entering  a  common  node 
must  be  zero.  Furthermore,  the  total  charge  is  assumed  to  be  distributed  over  each 
individual  wire  corresponding  to  its  wire  radius. 

For  the  ease  of  an  open-ended  wire,  one  cannot  expect  the  current  to  go  to  zero  at 
the  wire  end  if  the  wire  has  a  finite  radius.  For  a  ware  of  radius  a.  the  current  at  the 
wire  end  can  be  represented  as  a  function  of  its  derivative  as  [22 1 


/  ( z  )  I . 


-<f  ric  )  J  ,U'J  )  ji  (  -  ) 


:  al  end 
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where  u.  =  the  unit  normal  vector  of  the  end  cap  and  f  =  ±  1  indicates  the  direction 
ol  the  normal.  By  applying  the  current  and  charge  conditions  at  both  ends  of  each 
wme  segment,  two  of  the  three  coefficients  can  be  solved  for.  This  leaves  one 
i  oei  i  ic ic-nt  'o  be  solved  lor  by  the  matrix  inversion. 

I  he  .ncident  F-tield.  £'  .  on  the  left-hand  side  of  Equation  (2.1.2).  is  to  be 
d.c’a’ed  by  the  physical  problem.  In  this  applica’ion.  it  is  desirable  to  excite  a  thin- 
w,re  '.'Picture  with  a  vertically  directed  h  - 1  ic-  Id  distributed  over  the  surface  ot  a 
s;iui,!  region.  For  the  FMP  simulator  problem,  this  would  denote  the  source  gap  region. 
I  ne  \FC  otters  two  voltage  source  models.  1  he  tirst  is  the  Current  Slope 
Disc  on  t  .nuity  (  Si  > )  model,  w  hie  h  consists  ol  a  high  I  y  localized  voltage  source  across 
•  ne  bicm.ical  junction  of  two  wires,  [he  second  type  of  source  is  an  Applied  Electric 
(  c  ld  1  \FF)  source  model.  I  his  tinxiel  is  a  voltage  source,  of  constant  magnitude. 


(*sV 

>v. 


applied  over  a  designated  region  of  the  wire  structure.  The  resultant  electric  t  u  Id  >  ,m 
He  represented  as  a  pulse  t  unction,  whose  magnitude  can  be  expressed  as 


where  V',  is  the  complex  amplitude  ol  the  applied  voltage  on  'he  j  ~  scgnuni  and  A 
is  the  length  ot  the  / '"*  segment. 

1  he  AhF  is  a  desirable  source  model  for  defining  a  distributed  held  -uer  ”u 
source  gap  region  of  the  EMP  simulator  problem,  and  will  be  used  for  this  application. 
However,  one  must  realize  its  limitations  before  using  it  for  any  application.  The  Al  l 
is  constant  over  the  entire  wire  segment  it  is  applied  and  is  assumed  to  be  ero 
everywhere  else.  1  herefore.  on  the  left-hand  side  ol  the  matrix  lormulation  o| 
Equation  (2.1.2).  the  E-field  is  assumed  to  be  a  pulse  I  unction.  However  ’he  p 

discontinuity  of  the  field  is  nonphysical  and  the  actual  incident  held  m-.st 
calculated  from  the  computed  currents.  In  this  study,  it  is  found  that  under  cit'jip. 
conditions  the  field  can  vary  a  significant  amount  from  its  assumed  value. 
Understanding  this  is  quite  important  for  two  reasons.  I  he  lirst  is  the  desire  *o 
accurately  model  the  incident  E-field  source  for  a  reliable  approximation  ol  ”.e 
problem  being  solved.  Second,  the  input  impedance  of  the  source  is  calculated  bv  \E( 
assuming  a  constant  current  and  E-field  over  the  entire  source  region.  However  m 
i.er’a.n  applications,  this  may  not  be  a  good  approximation,  and  the  result  couid  ”)c 
misleading.  The  reader  is  relerred  to  the  paper  presented  by  Halpern  and  Mittra  (25 1 
Am  a  complete  analysis  on  the  validity  of  the  AEF  and  (A  D  voltage  source  nuxJels  tor 
'he  case  ol  ,i  thm-wire  antenna.  It  was  lound  that  lor  source  regions  much  smaller 
than  'he  leng'h  of  the  antenna,  and  tor  'he  radius  of  the  wire  sufficiently  small  l  <  < 
length  ol  the  source),  the  calculated  AEF  is  much  like  the  desired  pulse.  I  his  s  'rue 
since  the  thm-wire  approximation  being  used,  and  the  smaller  the  wire  radius,  'he 


better  the  approximation. 


2.2  Solution  ol  t he  Scattered  Fields  ol  I  hm-Wire  Structures 

Be!  ore  i"t p' i nji  the  BMP  simulator  problem,  a  much  simpler  example  w  ii  ■> 
analyzed  lirst.  It  a  ill  then  be  |>ossible  to  characterize  the  simulator  problem  b\ 
realizing  certain  limitations  encountered  when  using  NEC.  Therefore.  ’his  section  wi,; 
emphasize:  (  1  )  The  accuracy  ol  the  solution  due  to  the  number  ol  unknowns,  or  wire 
segments.  (2)  1  he  validity  ol  the  source  mixlel  due  to  the  characteristics  ol  trie 
antenna  model,  HI  The  determination  ol  conditions  necessary  to  render  the  solution 
invalid.  Three  examples  are  analyzed  lor  ’his  purpose:  l  I  )  Thin  wire  monopole  over  a 
PTC  ground  plane.  ( 2 )  Two-w  ire  rhombic  antenna,  t  $)  Planewave  scattering  ol  a  thin 
PEC  plate.  These  three  examples  are  presented  since  they  directly  relate  to  the 
accuracy  oT  carious  aspects  ol  the  parallel-plate  EMP  simulator  model.  I  he  thin-ware 
monopole  helps  >o  characterize  the  source  model  and  the  accuracy  of  the  solution  due 
>o  'he  number  ol  basis  t  unctions  chosen.  The  rhombic  antenna  problem  is  treated  as  a 
’ ransmission-line  type  ol  EMP  simulator  and  enables  one  to  evaluate  the  construction 

■  'I  ’he  source  and  load  regions  as  well  as  aid  in  the  interpretation  ol  results  ol  the 
parallel  plate  EMP  simulator.  F  inally  ’he  scattering  of  the  PEC  plate  enables  the 

■  haractenzation  ol  the  validity  of  the  solution  for  the  thin-wire  mesh  approximation 
'I  a  surt ace  when  subject  to  an  incident  t  ield. 

2  2  I  I  hin  wire  monopole 

I  v-  .'.*2  I  .11  .strates  the  'hin  wire  monopole  considered  in  this  section.  !  he  M  < 

,s  .sed  s« •  i v *•  i  or  'be  ax.al  i  ;rrer.»  produced  by  a  source  excitation.  Tne  axia. 

.rrent  .s  represented  by  the  'hree -term  nasis  t  unctions.  Recall  'hat  'he  U  s'iny 
'  unctions  ised  »o  genera’e  the  impedance  ma'rix  .it  'he  moment  •ne’hod  solut.on  we  re 
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chosen  to  be  mathematical  point  functions.  However,  when  using  point  matching  one 
has  to  try  to  avoid  sampling  error,  which  leads  to  an  incorrect  approximation  of  the 
current.  Therefore,  it  is  important  to  perlorin  a  boundary  condition  check  on  the 
a x id  1 1  >  directed  surface  E-f held  to  guarantee  that 

F!  =  -E!  (2.2.1) 

;s  satisfied.  It  is  recommended  that  at  least  10  samples  be  taken  per  wavelength  in 

order  to  minimize  sampling  error. 

First,  consider  the  case  ot  a  monopole  that  has  a  height  H-0.47  \  and  radius  a. 
where  H/a  -  1000.  Figure  2.2  illustrates  the  resultant  scattered  axial  E-field 
computed  along  the  antenna  surface  in  the  axial  direction.  Note  that  there  is  boundary 
condition  error  at  the  trailing  edge  of  the  source  region,  and  at  the  free  wire  end  of  the 
monopole.  The  singularity  due  to  the  free  wire  end  is  nonphysical,  and  is  a  result  of 
the  thm-wire  approximation. 

In  »he  study  of  most  antenna  problems,  it  is  desirable  to  compute  the  input 
mpedance  in  order  to  best  design  a  matching  network  for  the  antenna  excitation. 
From  the  reciprocity  theorem,  it  is  easy  to  develop  the  relationship  between  the 
'ermir.al  voltage  and  the  source  field  as  (2b) 

V.„  /„=/£'(.-)/(.’  )d:  (2.2.2) 

where 

*  =  the  JXldl  J  licet  H  >n 
F  =  \I  F  sou/cp 
/  (  :  )  =  [hr  j.xidl  i  ur  rent 

and  ’he  integration  is  performed  over  the  source  region.  Furthermore,  by  realizing 
'  hat  t  —  —h 

/  =  =  -JL  f  PI;  )H:  )J:  (2.2.1- 

l-  TJJ 

*huh  represents  'he  antenna  input  impedance  m  terms  of  the  computed  scattered 


It  can  be  shown  that  the  above  expression  is  stationary  and  is  very  insensitive 
to  small  errors  in  the  computed  current  (25).  As  the  limit  of  the  wire  segment  si/e 
becomes  ver\  small,  the  current  becomes  nearly  constant  on  the  wire.  I  herefore. 
assuming  l(  i )  -  a  constant  -  /  ,  .  the  above  integral  equation  reduces  'o 

1  in 

where  the  integral  is  over  the  source  region.  Now.  if  there  is  little  boundary  condition 
error  on  the  wire  sur! ace.  then  it  can  be  assumed  that 

=  -V',,  =  /  EH:  )J:  .  (  2.2.5  i 

f  here! ore,  using  the  expression  given  in  (2.1.6),  the  input  impedance  ol  the  antenna 

can  be  approximated  as 

V.„ 

An  =  y—  (  2.2.6  1 

1  n 

w  here 

V.,  =  in  pul  source  v<>  Itjge 

and 

l  „  -  the  current  ul  the  center  nf  ihe  wire  segment. 

!  his  x  -he  particular  approximation  used  by  NFC.  In  this  formulation  the  scattered 
held  -x  hich  is  computed  from  the  axial  wire  current,  is  assumed  to  be  a  pulse 
•  re’. on.  such  ’hat  the  field  is  constant  over  the  input  wire  segment,  and  zero 
-  v»-rv  w  nere  else.  Figure  2. 1  illustrates  the  act  ual  scattered  held  and  the  axial  c  urren’ 
i-  mpu'ed  over  'he  source  region  lor  the  example  considered  earlier  'Figure  2.2)  and 
pros  dm  a  comparison  with  ’he  AFF.  I  he  scattered  !  ie Id  deviates  somewhat  I  rom  the 
i  mxi  'ol  pulse  function,  which  Pads  to  a  small  error  in  the  input  impedance. 
However  tor  antennas  with  much  larger  radii,  ’he  scattered  field  deviates  greatly 
'  - ’m  ’he  pulse  *  um  'on  and  'he  approximate  expression  tor  'he  inpir  mi pedanc e  ,s  -o 


longer  valid.  In  this  case,  the  stationary  form  of  Equation  (2.2.3)  must  be  relied  upon. 
A  thorough  analysis  of  this  problem  is  provided  in  (25).  In  that  analysis  it  was  found 


that  if  the  radius  of  the  wire  was  on  the  order  of  one-thousandth  of  the  length  of  ihe 


monopole  antenna,  the  boundary  condition  error  became  very  small,  and  the 


approximate  form  of  the  input  impedance  was  accurate  to  within  a  small  percentage. 


However,  for  thicker  antennas,  where  the  radius  is  on  the  order  of  one-hundredth  of 


the  length  of  the  antenna,  the  approximate  form  is  not  as  accurate,  and  the  stationary 


form  of  the  input  impedance  should  be  used.  Furthermore,  for  the  thin-wire 


approximation  to  be  valid,  one  should  not  use  wire  segments  that  have  radii  larger 
than  one-eighth  of  the  segment  length:  otherwise,  the  approximation  breaks  down. 


There  is  one  important  limitation  encountered  when  using  NEC  that  is  not 


considered  by  (25).  and  that  is  the  problem  of  over  sampling.  In  other  words,  what 
happens  if  the  segment  size  becomes  too  small9  The  following  is  an  analysis  of  this 


problem. 


2.2. 1.1.  The  problem  of  over-sampling.  Over-sampling  can  occur  if  too  many  wire 
segments  are  chosen  to  represent  the  geometry  such  that  the  segment  size  becomes  too 


■jmail.  This  leads  to  unreliable  results  for  the  axial  current,  and.  hence,  the  scattered 


idds.  I  he  reason  for  this  error  is  purely  numerical.  As  discussed  earlier,  the  current 


s  expanded  into  a  series  of  three-term  current  basis  functions.  The  current  expression 
:s  mven  by  Equation  (2.1.4).  Initially.  A,  is  chosen  to  be  equal  to  -1  to  reduce  the 


iber  of  computations.  Then,  al  ter  satisfying  Kirchoff's  law  and  charge  continuity, 


’he  complete  basis  function  is  defined  by 


A,  =  -  I 


B  =(u,-g,- 


sin  *  A,  /  2 


sin  k  A, 


>.y  v'-x 


IX 


where 


C.  =  U~Q;"  -a:+Q,  +  ) 


cos  k  A,  /2 
sin  £  A, 


(22.7) 


u,  *(  1  —  <  <>s  k  A  )  —  P  *sm  k  A. 

<>  ~  =  _ : _ : _ : _ 

P  P  *  -eu(  cj  '  )  sin  A'  A,  +  (  P.  a,~  —  P  *al  )  cos  k  A. 

u.'Uos  k  A 1 )  —  P.  ~sin  k  A. 

o  *  = _ : _ ! _ : _ : _ 

(  P  P,~  -eu,  a,  *  )  sin  k  A,  +  ( Pt  at”  —  Pt  *u,  )  cos  k  A, 

and 


■'  1  —  cos  k  A , 

P.~  =  T  - L  up 

“T,  sm*  A, 

is  the  summation  of  the  V“  segments  connected  to  end  I  of  the  i‘h  segment. 

Furthermore. 


P:'=  I 


■  =1 


cos  k  A ,  —  1 


sin  k  A , 


is  the  summation  of  the  .V*  segments  connected  to  end  2  of  the  t‘n  wire  segment. 
Finally,  u,  contains  'he  unkown  coefficients,  and  is  a  function  of  the  segment  wire 
radius.  However,  it  is  important  to  understand  what  happens  to  the  <2, *  as  A, 

2  trA, 


becomes  much  smaller  than  A.  or  as 


becomes  very  small.  Analysing  the 


expressions  for  Pt  *  and  P  shows  that  they  area  function  of 


1  —  cos  k  A , 

±  _ _  . 

sin  k  A; 

Therefore,  as  A  —  0,  both  1  -cos  k  A,  0  and  sin  k  A .  -»  0.  Similarly,  in  Q,  *,  both 
The  numerator  and  denominator  tend  to  zero  as  the  segment  length  becomes 
increasingly  smaller.  Numerically,  a  computer  cannot  handle  this  and  the  program 
will  sutler  t  rom  underflow  during  the  matrix  fill.  Subsequently,  serious  errors  will 
result  in  the  computed  current.  It  should  be  noted  that  the  exact  place  that  underflow' 
will  occur  depends  on  the  precision  of  the  computer  used.  The  following  results  were 
obtained  for  NEC  written  in  double-precision  Micro- Soft  Fortran,  which  has  32  bytes 


of  precision,  and  running  on  an  AT&T  PC6300  desktop  computer. 

In  Volume  II.  Section  2  of  (22).  it  is  recommended  that  the  segment  length  should 
not  be  smaller  than  IO-JA  due  to  this  underf  low  error.  However,  the  program  itself 
does  not  l  lag  ’he  user  it  underflow  does  occur,  and  an  erroneous  solution  does  result. 
Therefore,  the  user  must  be  very  careful  in  interpreting  his/her  results.  The 
following  is  a  good  example  of  this  over-sampling  problem. 

Consider  a  short  monopole  over  a  PEC  ground  plane,  where  H/A  -0.1.  and  H/a  = 
1000.  Table  2.1  illustrates  the  self-term  of  the  matrix,  computed  over  the  source 
segment,  as  well  as  the  current  and  surface  E-field,  as  a  function  of  the  number  of 
wire  segments.  Note  that  the  AEF  is  applied  over  the  first  segment  only  in  each  case, 
and  all  segments  are  of  equal  lengths.  For  this  example,  it  was  found  that  if  the 
segment  length  becomes  <  0.005  A,  the  solution  becomes  quite  unreliable.  Therefore, 
when  the  segment  length  is  on  the  order  of  10"J  A.  it  is  very  important  to  check  the 
accuracy  of  the  solution  with  at  least  a  boundary  condition  check  of  the  axial  surface 
E-f ield.  or  by  monitoring  the  self -terms  of  the  short  wire  segments. 

2. 2. 1.2.  The  source  model.  This  thesis  is  concerned  with  the  time  domain  pulse 
excitation  of  an  antenna.  Therefore,  when  solving  for  the  currents  and  scattered  fields 
in  the  frequency  domain,  it  is  necessary  to  solve  the  problem  over  a  broad  frequency 
range.  However,  for  each  frequency  it  is  desirable  to  minimize  the  number  of 
unknowns  to  reduce  the  computation  time,  but  still  obtain  an  accurate  solution.  In 
order  to  model  the  antenna  in  a  consistent  manner,  the  source  region  must  be  the  same 
tor  all  frequencies.  Therefore,  the  segments  denoting  the  source  region  must  be  small 
enough  such  that  the  wire  segments  do  not  exceed  0.1  A  in  length.  This  can  easily  be 
avoided  by  cascading  voltage  sources  in  series  if  necessary.  However,  one  must  also  be 
caret ul  to  ensure  that  the  source  does  not  become  too  small;  otherwise,  the  result  will 
suffer  from  over-sampling  error.  This  could  be  a  severe  limitation  when  considering 


problems  where  the  wavelength  is  much  longer  than  the  source  region.  For  much 
larger  problems,  such  as  the  parallel-plate  bounded-wave  simulator,  this  obviously  is 
much  more  of  a  concern  since  the  source  region  is  much  smaller  than  the  size  of  the 
■or  ict  ire  Finally,  one  must  be  careful  ’hat  the  wire  radius  does  not  exceed  one-eighih 
o!  'he  wire  length;  otherwise,  the  thm-wire  approximation  will  break  down. 

With  the  use  of  the  monopole  antenna  problem  some  intuitive  insight  to  the 
n’erpretation  ot  the  results  of  the  larger  HMP  simulator  problem  can  be  gained.  T  he 
rases  of  the  antenna  being  2  A  and  I  \  in  length  are  considered.  In  the  example,  the 
source  region  will  be  composed  of  2  wire  segments,  where  the  AEF  is  distributed  only 
on  the  first  segment,  and  is  of  height  HS.  The  effects  on  the  segment  current  as  the 
source  region  becomes  very  small  is  illustrated  in  Table  2.2.  As  was  found  previously, 
error  in  the  self-term  arises  when  the  length  of  the  segment  is  on  the  order  10-3. 
Figure  2.4  illustrates  the  axially  directed  E-field  on  the  wire  surface  over  the  source 
region  as  HS  reduces  in  size.  A  significant  boundary  condition  error  is  noted  m  Figure 
2.4(c).  when  HS  IK  -  self-term. 

2.2.2  Rhombic  EMP  simulator 

The  parallel-plate  bounded-wave  EMP  simulator  is  usually  referred  to  as  a 
transmission  line  type  of  simulator  since  it  supports  guided  waves,  and  is  typically 
•errninated  by  a  matching  impedance.  Therefore,  the  fields  can  be  expressed  as  a 
viperposition  of  discrete  and  continuous  modes.  V  arious  modal  solutions  for  the  finite 
■j.  :d  t  h  plate  geometry  are  discussed  in  Chapter  1.  However,  these  analogies  do  not 
'real  'be  problem  as  a  finite  length  transmission  line,  but  as  an  infinite,  or  semi- 
.niimte  proniem.  In  order  to  develop  a  better  understanding  of  the  parallel-plate 
bounded-wave  EMP  simulator  as  a  transmission  line  type  of  simulator,  it  will  be 
instructive  ’o  analyze  an  elementary  configuration  first.  In  order  to  do  this  the 


physical  attributes  of  the  simulator  geometry  is  utilized. 

As  a  result  of  the  finite  width  of  the  conducting  plates  of  the  parallel-plate 
simulator  most  of  the  current  will  be  distributed  at  the  edges  ol  the  plate.  Therefore, 
j  l  .rst-'rder  approximation  would  be  to  limit  the  current  to  a  line  current  located  at 
the  outer  edges  of  the  plates.  Furthermore,  if  the  line  currents  follow  the  edges  of  the 
conical  end  sections  and  meet  in  the  center  above  the  parallel-plate  region,  the 
geometry  resembles  that  of  the  well-known  rhombic  antenna.  The  far-lield 
characteristics  of  the  two-wire  conductor  rhombic  antenna  are  well  known  [27); 
however,  it  is  the  near-field  characteristics,  computed  between  the  conducting  wires, 
that  are  of  interest  for  the  EMP  simulator  application.  More  recently,  Shen  and  King 
have  identified  the  usefulness  of  the  rhombic  antenna  as  a  transmission-line  type  of 
EMP  simulator  (28-32).  Not  only  did  they  show  the  fields  to  be  quite  uniform 
between  the  conducting  wires,  and  predominately  TEM,  bui  they  also  identified  its 
usefulness  as  an  analogous  problem  to  the  parallel-plate  type  of  simulator.  In  their 
works,  an  extensive  experimental  and  theoretical  analysis  of  the  problem  is  provided, 
and  the  reader  is  encouraged  to  review  them. 

2.2  2.1.  Numerical  solution  using  NEC.  In  this  section,  the  solution  of  the  field 
distribution  within  the  bounded  region  of  the  conducting  wires,  when  the  rhombic 
EMP  simulator  is  excited  by  a  continuous-wave  excitation  across  the  source-gap  region 
and  'erminated  by  a  purely  resistive  load  impedance,  is  derived  numerically  using 
NEC  The  results  are  then  compared  to  the  theoretical  and  experimental  results 
published  by  Shen  and  King  [28- 12). 

I  he  rhombic  antenna  is  modeled  as  a  thin-wire  structure,  and  the  currents  on  the 
w  ires  are  assumed  by  the  three-term  basis  function  described  earlier  in  this  chapter. 
However,  complications  arise  while  trying  to  establish  a  proper  source  and  load  model 
tor  The  t ransmission-line  type  of  simulator.  Therefore,  an  analysis  of  this  problem 


will  be  extremely  useful  in  establishing  a  model  of  the  source  and  load  regions  that 
successfully  approximates  previous  theoretical  and  experimental  results.  These 
results  will  be  used  in  the  analysis  of  the  parallel-plate  simulator  problem.  The 
following  is  the  development  of  the  source  and  load  models  using  a  thm-wire 
approximation. 

Source  Model. 

The  source  model  is  extremely  important  in  establishing  the  proper  excitation  of 
the  transmission-line  type  of  EMP  simulator.  The  fields  produced  in  the  region 
between  the  conducting  wires,  usually  referred  to  as  the  working  volume,  are  expected 
to  be  TEM  for  lower  frequencies  as  the  structure  is  essentially  a  simple  transmission 
line.  The  source  must  be  oriented  such  that  the  TEM  mode,  which  contains  a  vertical 
E-field  component,  is  excited.  Therefore,  it  is  desired  source  is  a  time  harmonic 
potential  of  constant  amplitude  established  by  a  vertical  E-field  distributed  across  the 
2-dimensional  source  gap.  However,  using  a  thin-wire  antenna,  modeled  by  NEC.  one 
is  limited  to  a  finite  length  line  source.  The  Applied  E-Field  ( AEF)  source  discussed  in 
Section  2.1.  which  is  a  constant  E-field  distributed  over  a  wire  segment,  best 
approximates  the  desired  excitation.  Figure  2.5(a)  illustrates  the  thin-wire  model  of 
the  source  region.  The  applied  voltage  is  expressed  as 

v  =Je‘-jT. 

It  must  be  realized  that  there  is  a  definite  tradeoff  to  this  configuration  due  to  the 
introduction  of  a  sharp  bend  in  the  thin-wire  geometry.  In  Section  2.1  it  was 
discussed  that  charge  and  current  continuities  are  established  at  the  junction  of  two  or 
more  wires.  However,  even  though  the  continuity  of  total  current  is  maintained,  the 
individual  basis  vectors  of  the  current  are  discontinuous.  This  discontinuity  in  the 
current  leads  to  a  nonphysical  singular  behavior  in  the  E-field,  resulting  in  a 
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boundary-condition  error  of  the  tangential  E-field  on  the  wire  surface.  Therefore,  one 
must  be  careful  that  this  singularity  does  not  significantly  affect  computed  current 
distribution  on  the  conducting  wires  throughout  the  structure  or  the  characteristics  ot 
the  simulator. 

Figure  2.6  illustrates  the  computed  scattered  E-field  along  the  axial  direction  on 
the  surface  of  the  source-gap  wire  segment  due  to  a  constant  incident  E-tield 
distributed  over  the  length  of  the  wire.  The  boundary  condition  error  is  apparent  at 
the  bend:  however,  it  was  found  that  this  configuration  yielded  the  best  results  in 
comparison  to  Shen  and  King's  data. 


Terminating  Load  Impedance. 


In  order  to  avoid  standing  waves  on  the  transmission  line,  the  rhombic  simulator  j 

must  be  terminated  in  its  characteristic  impedance.  For  low  frequencies,  where  the 
center  height  is  much  smaller  than  a  wavelength,  the  characteristic  impedance  is 
nearly  constant  and  the  simulator  is  best  matched  by  a  purely  resistive  load.  The 
dominant  mode  is  TEM  to  y.  and  the  E-field  is  in  the  vertical  direction,  z.  Therefore. 

'he  load  should  be  oriented  in  the  same  direction  as  the  vertical  E-field  such  that  the  ] 

energy  is  dissipated  by  the  load.  Since  £  is  related  to  the  electric  vector  potential  as  , 

£  =  —  j  ai/ju\  +  VV'rA  I 

the  current  induced  on  the  wire  will  be  in  the  direction  of  the  incident  E-field. 

Therefore,  the  wire  load  should  be  oriented  in  the  z-direction.  The  geometrical 

coni  iguration  of  the  wire-load  is  similar  to  that  of  the  source,  where  the  loading 

resistance  is  distributed  uniformly  over  the  wire.  It  was  found  that  the  transmission  J 

line  was  best  matched  by  placing  the  load  in  parallel  across  the  load  gap.  This  is 

illustrated  in  Figure  2.5(b)  (page  45).  Again,  one  must  be  aware  of  the  existence  ol  > 

boundary-condition  error  due  to  the  abrupt  bend  in  the  wire. 
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Characteristic  Impedance. 


Il  is  expected  that  if  the  t ransmission-lme-t y pe  simulator  is  indeed  terminated  in 
its  characteristic  impedance,  then  the  current  along  the  wires  will  be  purely  'raveling 
wave.  However,  because  of  the  existence  of  higher-order  modes,  there  will  be  some 
reflection  from  the  load.  Furthermore,  the  characteristic  impedance  of  the  line  will 
change  with  frequency,  and  in  order  to  determine  its  value,  the  analogy  used  by  Shen 
and  King  [28]  will  be  followed.  The  Standing  Wave  Ratio  (SWR)  of  the  current  is 
found  as  a  function  of  the  load  impedance.  When  the  SWR  is  minimized,  the  line  can 
be  assumed  to  he  matched.  At  lower  frequencies,  the  simulator  fields  are  bounded  and 
most  of  the  energy  is  dissipated  by  the  load.  However,  small  deviations  in  load 
impedance  will  lead  to  large  changes  in  the  SWR  of  the  line.  At  higher  frequencies, 
the  structure  will  radiate  a  significant  amount  of  energy.  Therefore,  there  is  less 
energy  to  be  dissipated  by  the  load,  and  the  resultant  reflected  waves  will  be  smaller. 
This  is  illustrated  by  Figure  2.7.  At  lower  frequencies,  more  energy  is  distributed  in 
the  vertical  E-field,  which  is  predominately  TEM.  and  is  dissipated  by  the  load. 
However,  at  higher  frequencies,  higher-order  modes  exist,  and  a  larger  percentage  of 
the  energy  is  distributed  in  the  transverse  and  longitudinal  components  of  the  field, 
which  are  not  dissipated  by  the  load.  Therefore,  the  load  is  often  sloped  or  distributed 
such  that  these  higher-order  modes  can  be  dissipated,  and  not  reflected  back  into  the 
center  parallel-plate  region. 

2. 2. 2. 2.  Characteristics  of  the  rhombic  EMP  simulator.  The  purpose  of  this  analysis 
is  to  investigate  the  computed  characteristics  of  the  transmission-line-type  rhombic 
EMP  simulator  using  the  Numerical  Electromagnetics  Code.  Initially,  the 
characteristic  impedance  is  determined  by  finding  the  terminating  load  impedance  that 
minimizes  the  standing- wave  ratio  of  the  conducting  wire  currents.  The  currents  and 
the  input  impedances  are  analyzed  as  a  function  of  the  load  as  well  as  a  function  ol 
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frequency.  Finally,  the  electric  field  distribution  within  the  area  between  the 
conducting  wires  is  examined,  and  its  uniformity  as  a  function  of  frequency  is 
determined.  The  results  are  compared  to  both  the  experimental  and  theoretical  results 
published  b\  Shen  and  King.  There! ore.  the  same  simulator  geometry  studied  by  Shen 
and  King  [28]  is  used  here.  This  is  illustrated  in  Figure  2.8. 

The  broad-band  characteristics  ot  the  simulator  can  be  defined  by  three  regions  of 
interest.  In  the  low-frequency  region,  where  the  wavelength  is  much  larger  than  the 
maximum  height  of  the  simulator,  the  field  has  predominately  TEM  characteristics 
and  is  quite  uniform.  In  the  middle-frequency  region,  the  height  of  the  simulator  is 
on  the  order  ot  a  wavelength  and  higher-order  modes  become  more  significant. 
Finally,  in  the  high-l  requency  region,  the  structure  acts  as  a  radiator,  and  very  little 
energy  is  actually  dissipated  by  the  terminating  load  impedance.  For  simplicity,  each 
ot  the  three  regions  is  defined  by  a  single  frequency.  The  regions  are  defined  by:  1) 
x=4m.2)\=lm.  1 )  \  =  0.5  m  . 

Initially,  the  characteristic  impedance  of  the  rhombic  simulator  in  each  of  'he 
'hree  frequency  regions  is  determined.  Table  2.3  illustrates  the  Standing  Wave  Ratio 
i  SWR  )  of  the  current  as  a  function  of  load  impedance,  computed  from  the  NEC  model. 
Furthermore,  it  also  lists  the  efficiency  of  the  simulator,  which  is  a  measure  of  the 
percentage  ot  power  radiated  to  the  total  input  power  available.  When  the  line  is 
properly  matched,  the  efficiency  is  near  minimum  since  a  maximum  amount  of  power 
s  being  dissipated  by  the  load.  Therefore,  the  efficiency  can  also  be  used  as  a  means  to 
determine  the  characteristic  impedance.  Table  2.4  provides  a  comparison  of  the  results 
obtained  !  rum  ’he  NFC  approximation  with  those  obtained  by  Shen  and  King.  In  the 
middle-  and  !ow -I  requency  regions,  the  results  compare  fairly  well.  The  deviation  is 
mainly  due  to  the  source  and  load  model  deficiencies.  However,  in  the  high-frequency 
region  the  results  deviate  by  a  signif icant  amount.  This  is  due  to  limitations  associated 


with  the  source  and  load  models.  As  the  wavelength  becomes  comparable  with  the 
dimension  of  the  source  and  load  gaps,  the  approximation  seems  to  break  down. 
Although  it  is  interesting  that  even  though  the  matching  impedance  is  much  lower  m 
’he  high -t  req uenc v  region,  the  near  holds  and  the  conducting  wire  currents  still 
compare  quite  well  with  those  obtained  by  Shen  and  King. 

It  should  be  noted  that  even  though  the  characteristic  impedance  changes  wor. 
frequency,  the  SWR  of  the  mismatch  throughout  the  band  ot  interest  can  be  less  than 
1.2.5  when  matched  to  the  low  frequency  characteristic  impedance.  Furthermore  o 
should  be  realized  that  under  pulse  excitation,  most  of  the  energy  typically  lies  in  the 
low-  to  mid-frequency  region.  This  is  encouraging  since  it  means  that  a  load 
impedance  can  be  developed  that  will  efficiently  terminate  the  simulator. 

Figure  2.9  provides  a  comparison  ot  the  wire  currents  ot  the  open-circuited 
simulator  computed  by  NEC.  with  those  published  by  Shen  and  King,  when  A.  =  4  'n 
Figure  2.10  compares  the  computed  wire  currents  for  the  short  circuit  case,  when 
a  =  1  ”1  .  In  each  case,  the  NEC  solution  compares  quite  well  with  Shen  and  King's 
results.  However,  a  slight  shift  in  the  peaks  of  currents  provided  by  the  NEC  model 
-mould  be  noted.  The  amount  of  shift  is  equal  to  the  height  of  the  source  and  load 
-egions.  Therefore,  in  the  high-frequency  region,  when  the  gap  height  becomes  more 
sigml  icant  in  terms  of  a  wavelength,  large  errors  can  occur  in  the  computed  matched 
terminating  impedance. 

Figures  2.11  -  2.13  illustrate  'he  currents  along  the  conducting  wire  axis  as 
computed  by  NEC.  These  results  can  be  compared  to  Shen  and  King's  work 
. !  1  r.i'ed  m  f  igure  2.14.  and  are  found  to  compare  quite  well  m  all  frequence' 
regions. 

f  igures  2  15  2.2  3  illustrate  'he  F  fields  in  the  center  region  ol  the  simulator 

•  imputed  tor  matched  terminations  by  both  the  NEC  model,  and  by  Shen  and  King.  It 
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should  be  noted  that  in  the  experimental  results  presented  from  Shen  and  King's  work 
[2d],  when  the  geometry  is  scaied  down  by  i/4,  the  operating  wavelengths  are  staled 
dow  n  as  well.  Again  the  results  are  m  good  agreement. 

At  ios  t  rcquc-r.cies.  the  magnitude  of  the  vertical  E-t  ie!d  is  quite  uniform  w  /n.n 
•he  center  region  and  has  an  SWR  very  close  to  one.  Furthermore.  F .  is  a  progressive 
wave  with  a  linear  phase  shift.  In  the  middle-  and  high-l  requency  regions,  'he 
vertical  I, elds  are  again  quite  uniform  and  continue  to  maintain  a  nearly  linear  phase 
sh.tt  over  distance.  However  due  to  the  existence  of  higher-order  modes,  the  SW  k  oi 
the  tield  is  much  larger,  and  is  approximately  1.25  for  the  middle-f requency  region 
and  1.75  for  the  high-t  requencv  region.  It  should  also  be  noted  that  the  longitudinal 
component  of  the  E -field.  fv  .  is  of  significant  amplitude,  although  it  is  still  less  than 
the  veriical  component. 

In  conclusion,  good  agreement  was  found  between  the  Numerical 
Electromagnetics  Code  s  model  of  the  rhombic  EMP  simulator  and  the  theoretical  and 
experimental  results  published  by  Shen  and  King.  The  source  and  load  models  were 
developed  using  a  thm-wire  model  and  a  distributed  E-field  source  and  resistive-load, 
respectively.  The  computed  characteristic  impedance  of  the  simulator  in  the  low-  and 
middled  requency  regions  compared  quite  well  with  experimental  data.  However  u  it 
*o  'he  wire  model  of  the  source  and  load  gap,  the  matching  load  impedance  suffered 
significant  error  in  the  high-frequency  region.  Finally,  in  all  frequency  regions  'he 
computed  currents  and  E-f iclds  compared  extremely  well  with  those  ot  Shen  and  Kmi? 

under  matched  conditions,  regardless  of  the  error  in  the  value  of  the  characteristic 
i  m  pedances. 


2.2.  f  Scattering  of  a  C  plate 


[he  c  on  cl  ui  ’  in  a  surl  ace  o|  'he  para!  lei -plate  simulator  is  to  ’*•  appro'  ima’ed  b\  a 
thin -wire  mesh.  Imtiallv  ,  ,t  will  be  mst  r  active  to  develop  a  similar  problem  sun  tha' 
the  accuracy  of  the  wire  mesh  approximation  and  the  limitations  assiKiated  with  ,t  can 
easily  be  studied.  However,  in  order  to  do  this  effectively,  the  physical  attributes  oi 
the  currents  induced  on  the  conducting  plates  of  the  simulator  must  be  unde.stcxxi. 
The  geometry  of  the  simulator  is  an  open,  f mite- width  transmission  line.  Due  to  Tu 
finite  width,  there  will  be  a  singularity  in  the  current  at  the  plate  edge.  Furthermore 
»he  current  will  be  highly  oriented  m  the  longitudinal  direction,  and  very  little 
current  will  actually  flow  in  the  transverse  direction.  In  this  section,  the  currents 
induced  on  a  PEC  plate  scatterer.  when  illuminated  by  an  incident  plane  wave  a:e 
studied  A  comparison  will  be  made  between  the  solutions  derived  by  NEC.  using  a 
wire  mesh  approximation,  and  that  derived  by  a  more  accurate  technique  that  solves 
tor  the  two-dimensional  surface  current  distribution.  The  PEC  plate  scatterer  has  a 
t  mite  w udth  and  the  dimension  will  be  chosen  to  be  on  the  order  of  one  wavelength  in 
order  to  reduce  the  number  of  unknowns.  Furthermore,  in  order  to  reproduce  some  ot 
■he  limitations  associated  with  the  parallel-plate  simulator’s  wire  mesh  approximation 
r fie*  induced  currents  are  predominately  oriented  in  one  direction.  Therefore,  the  case 
ot  a  norma  iiy  incident  FE  plane  wave  excitation  is  considered. 

2.2.  t .  1  Numerical  Results.  Figure  2.21  illustrates  the  geometry  of  the  PEC  plate.  I  he 
plane  wave  incident  is  IE  to  y.  with  the  E-field  directed  along  the  /-direction.  1  he 
magnitude  ot  the  tangential  H-tield.  //,  ,  is  chosen  to  be  unity.  Therefore.  I  rom  a 
geometrical  optics  point  of  view,  the  induced  elect ric  currents  are  expected  to  be  on  t he 
order  ot  one  for  a  normally  incident  excitation,  f  he  plate  is  chosen  to  be  square  with 


a  width  of  1.1  wavelengths. 


I  —  '  ally  on  >1  'he  ’hm  I’f  (  plate  problem  is  obtained  using  a  xpec'ral 

(ia;erk.n  nr<x  ielurt‘  lommned  won  'be  me  t  bixl  ot  (  onjugjte  (gradients.  I  he 
T'  h  eo.  . '  A  :s  .oped  ->v  Rav  and  Mittra  a’  'be  I  mversitv  ol  Illinois,  and  is 

Uese'OVd  .n  it 'an  n  I  O)  (be  '  w  o  •dimensional  current  is  expanded  in  ’errns  .>1  a 
general  set  oi  nasis  :  unctions,  ret  erred  *o  as  triangular  patch  basis  l  unctions.  These 
"msis  tunc’ions  were  tost  considered  by  Clisson  in  [33).  The  basis  Junctions  are 
cien trails  spaied  such  'hat  the  FKI  operator  can  be  used  to  compute  'be  discretized 
..Vegral  operator.  I  hev  also  guaran'ee  'hat  the  current  normal  to  an  interior  edge  ot 
'he  patch  is  continuous  and  does  not  generate  any  nonphysical  line  cnarges. 
f  urthermore  only  parallel  components  ot  the  current  will  be  nonzero  on  the  edge  of 
'he  plate.  The  1.1  wavelength  pla'e  was  discretized  into  an  array  of  lb  by  lb  patches, 
each  containing  three  basis  (unctions,  leading  to  7b8  unknowns.  Convergence  to  less 
'han  0.1*)  error  required  only  b8  iterations.  However,  since  the  source  is  highly 
oriented  in  'he  z -direction,  the  current  ,n  the  x-direction  still  had  considerable  error, 
f  g  ares  2  22  and  2.2  f  illustrate  the  surtace  current  densities  in  'he  z-  and  x-directions. 
'especWeiv  a( 'er  TOO  iterations.  As  expected,  the  parallel  current  along  the  edges  ol 
'he  patch  is  singular.  Fur' hermore  7,  which  .s  'he  current  perpendicular  to  the 
;xiiari.'at ic.m  ol  'ne  incident  E-field.  is  approximately  three  orders  c>t  magnitude  less 
'  nan  7.  . 

•x  : * h  the  use  ot  NEC.  the  patch  surtace  is  approximated  by  a  wire  mesh.  As  was 
d.scussed  irrSection  2.1.  at  least  10  basis  functions  per  wavelength  must  be  used  ,n 
■  ret' r  ’»)  t-nsure  'he  proper  convergence  ol  the  solution.  However  it  should  tirst  be 
realized  that  'here  are  limitations  associated  with  'he  wire  grid  representation. 
Im’idil v.  'he  currents  are  limited  to  a  one-dimensional  flow  along  individual  wares, 
and  are  coupled  to  transverse  wires  at  multiple  wire  junctions.  In  Sec’ion  2.1  it  was 
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discussed  that  at  each  wire  junction  Kirchoff's  current  law  and  charge  continuity  are 
oemg  applied  as  boundary  conditions.  However  This  results  m  a  step  discontinuity  in 
'he  current  at  catn  junction.  Since  /.  .s  the  dominant  current,  its  discontinuity  is 
■elativiic  '.'Mall.  However  Jt  is  much  smaller  in  magnitude.  As  a  result,  the 
discont mud  v  in  the  current  is  quite  large. 

The  1.1  wavelength  patch  surface  was  approximated  by  a  12  by  12  square  grid, 
eadmg  'o  a  total  of  288  unknowns.  Figures  2.24  and  2.25  illustrate  the  surtace 
..urrent  densd.es  in  the  z-  and  x -directions.  respectively.  We  find  that  J.  compares 
extremely  well  in  magnitude  with  the  PEC  plate  problem.  The  singularity  in  the 
current  along  the  plate  edge  is  apparent;  however  it  is  not  nearly  as  significant  as  that 
computed  by  the  iterative  procedure.  Furthermore,  there  is  significant  error  in  the  x 
component  ol  the  current  and  this  resuit  is  unreliable,  although  it  cannot  be  thrown 
away  since  d  does  contribute  to  J .  at  the  wire  junctions.  Figure  2.26  illustrates  tr.e 
comparison  of  the  NFC  and  iterative  solutions  of  J.  along  the  z-direction.  Near  the 
plate  edge  the  normal  current  should  vanish.  As  discussed  earlier,  this  is  ensured  by 
•he  choice  ol  hasis  I  unctions  in  the  iterative  procedure;  however,  in  the  NEC  solution 
'he  wire  segment  terminates  into  a  multiple  wire  junction,  and  is  nonzero  because  d 
.st  satisfy  Kirchoff's  law  at  this  point. 

l  .giire  2  26  illustrates  the  current  distribution  ol  the  NEC  solution  located  at  tn.e 
,-n'er  'I  eacn  wire  segment.  However  Figure  2.27  illustrates  the  current  computed 
...mg  ’ne  length  ol  each  segment.  I  his  c  lear'v  shows  the  discontinuity  that  exists  m 
ict'ent  which  is  a  result  ol  the  mul’iple  wire  junctions.  Recall  that  thecurren’ 
n  •  ne  >  o.r e<  Mon  s  on  the  order  ol  a  magnitude  smaller  than  the  z  directed  current. 
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J,  s  quite  significant.  I  his  is  illustrated  in  Figures  2.28  and  2.2'* 
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Furthermore,  the  discontinuities  in  the  currents  lead  to  nonphysical  singularities  in 
'he  F- tie  Id  near  the  wire  mesh  surface. 

F he  analysis  of  the  parallel-plate  HMP  simulator  will  require  the  computation  ot 
'  he  near  - ;  .t  !d  u  ,st  r:  hut  ion  m  the  volume  directly  underneath  t  he  top  conducting  plu'c. 
Iherefore.  wnen  'he  conducting  surface  is  approximated  by  a  wire  mesh,  the  amour, ' 
.'t  distortion  of  the  H-field.  due  to  these  singularities,  must  be  determined.  Therefore 
n  C  hapter  ^  an  analysis  of  smoothing  of  the  current  along  both  directions  will  be 
nert  ormed.  and  the  ef  t  ect  on  'he  near  fie  ids  will  be  presented. 


TABLE  2.1  OVERSAMPLING 
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TABLE  2.2  ANTENNA  CHARACTERISTICS  DUE  TO  AEF  SOURCE  HEIGHT 
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TABLE  2.3(a) 

NEC  MODEL  OF  THE  RHOMBIC  EMP  SIMULATOR:  \  =  4  m 


z, 

Re(Z,„  ) 

Im(Z„  ) 

S1V7? 

(  Q  ) 

(O  ) 

(  O  ) 

(%) 

0.0 

38.80 

-20.28 

100.0 

1  1.60 
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73.22 

-22.54 

40.28 

4.20 

100.0 

104.80 
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2.56 
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-34.26 
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mo 

203.80 
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268.20 
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32.64 

2.46 

10000.0 
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TABLE  2.3(b) 

NEC  MODEL  OF  THE  RHOMBIC  EMP  SIMULATOR:  A  =  I  m 
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42.20 
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166.2 
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1.13 
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TABLE  2.4  CHARACTERISTIC  IMPEDANCE  OF  THE  RHOMBIC  EMP  SIMULATOR 
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Current  Along  Wire  Axis  (74.95  MHz) 


'omputed  by  Shcn  and  King  (29j. 


(a)  Magnitude  of  axial  current. 

Normalized  current  computed  along  the  conducting  wires  of  the  rhombic  fcMP 
imulator  by  NHC  under  matched  conditions  when  A  =  4  m  . 
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(a)  Virtual  I  I  it' Id  magnitude  produced  by  a  unit  Al:l;  excitation,  Computed  along  the 
ground  plane  in  the  longitudinal  direction,  y.  under  matched  load  conditions  for 
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ligure  2.17.  ( b)  Vertical  I:  ( leld  phase  couiputecl  along  the  ground  plane  in  the  longitudinal  direction 


Vertical  h-field  magnitude  computed  at  various  heights  in  the  working  volume  in 
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fb)  Iheoretical  and  experimental  vertical  E-fields  on  the  ground  plane 
computed  along  the  longitudinal  direction  [29).  Dimensions  are  scaled 
down  by  3/4. 

Figure  2.18.  Theoretical  and  experimental  results  of  the  distributed  E-field  in  the  x-0 
plane,  presented  by  Shen  and  King  for  K  =  1  m  . 
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Vertical  I:  field  magnitude  produced  by  a  unit  Ahl-  excitation.  Computed  along  the 
und  plane  in  the  longitudinal  direction,  y.  under  matched  load  conditions. 
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(a)  Theoretical  vertical  E-field  computed  along  the  longitudinal  direction 
[28]. 


(b)  Theoretical  and  experimental  vertical  E-fields  on  the  ground  plane 
computed  along  the  longitudinal  direction  [29].  Dimensions  are  scaled 
down  by  J/4. 


Figure  2.20.  Theoretical  and  experimental  results  of  the  distributed  E-field  in  the  x-0 
plane,  presented  by  Shen  and  King  for  A.  =  0.5  m  . 
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Figure  2.20.  (c)  Theoretical  longitudinal  E-field  computed  along  the  vertical  direction 
at  various  longitudinal  positions  [28). 
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CHAPTER  3 


NUMERICAL  ANALYSIS  OF  THE  PARALLEL-PLATE 
BOUNDED- W  AVE  EMP  SIMULATOR 

In  Chapter  l,  a  number  of  analogies  for  approximating  the  field  distributions 
within  the  parallel-plate  region  of  the  parallel-plate  bounded  wave  EMP  simulator 
were  presented.  These  approximations  are  all  limited  to  the  case  of  a  parallel-plate 
transmission  line  that  is  either  semi-inf intte  or  infinite  in  length.  Furthermore,  the 
approximations  are  typically  only  valid  for  low  frequencies,  and  for  the  strict 
parallel-plate  geometry.  However,  in  practice,  the  simulator  consists  of  a  finite  length 
parallel-plate  transmission  line  and  is  generally  excited  and  terminated  by  tapered  end 
sections.  In  this  chapter,  an  alternate  solution  is  considered.  The  simulator  is  treated 
as  a  large  conducting  body  and  the  finite  length  simulator  problem  is  solved  with  the 
use  of  scattering  theory.  This  approach  will  not  only  be  applied  to  the  typical 
parallel-plate  simulator  configuration,  but  will  also  be  applied  to  the  analysis  of  the 
simulator-obstacle  interaction. 

3. 1  Development  of  the  Scattering  Problem 

The  scattered  field  is  the  field  produced  by  the  induced  currents  and  charge 
distributed  on  some  body  when  illuminated  by  an  incident  field.  For  the  finite  length 
simulator  problem,  the  incident  field  is  taken  to  be  the  E-field  distributed  over  the 
source  gap.  The  scattered  field  is  then  found  by  satisfying  the  boundary  conditions  ot 
the  t-field  on  the  conducting  surface  of  the  simulator,  and  over  the  distributed 
vrmmatmg  impedance,  and  subsequently  solving  the  EFIE. 

As  was  discussed  in  Chapter  2.  the  EFIE  cannot  be  solved  analytically  in  general, 
and  a  numerical  solution  for  the  current  on  the  conducting  surface  must  be  developed. 
To  solve  the  integral  equation,  one  must  discretize  the  two-dimensional  surface 
current  into  a  truncated  series  of  known  basis  functions,  and  subsequently  solve  for 


Ihe  unknown  coefficients  by  moment  methods  or  an  iterative  procedure.  Although 
computer  memory  and  speed  are  becoming  more  accommodating  for  larger  problems, 
the  amount  of  programming  required  to  derive  a  solution  is  inexpedient.  Furthermore, 
the  size  of  the  simulator  is  generally  much  larger  than  the  wavelength  of  interest,  and 
the  number  of  unknowns  required  to  converge  would  be  too  large  to  justify  this 
approach. 

In  Chapter  2.  three  problems  leading  to  the  development  of  a  thin-wire 
approximation  of  the  conducting  surface  of  the  parallel-plate  simulator  were 
illustrated.  From  these  examples,  it  is  seen  that  NEC  can  be  a  very  useful  tool  for 
analyzing  the  parallel-plate  EMP  simulator.  Not  only  can  it  provide  a  good 
approximation  of  the  problem,  but  it  provides  an  easy  alternative  to  analyzing  the 
problem. 

The  conducting  surface  of  the  parallel-plate  simulator  is  two  dimensional  and  can 
easily  be  approximated  by  a  rectangular  wire  mesh.  This  configuration  ensures  both 
longitudinal  and  transverse  current  flows  along  the  conducting  surface.  The  size  of 
the  mesh  is  very  important  to  ensure  the  convergence  of  the  solution  and  should  be 
made  much  smaller  than  a  wavelength  in  order  to  better  approximate  the  surface 
currents.  However,  if  it  is  too  small,  it  will  increase  the  number  of  unknowns  and  the 
problem  becomes  too  large.  In  Chapter  2.  it  was  observed  that  at  least  10  basis 
functions  per  wavelength  were  required  in  order  to  properly  converge  to  a  solution. 
Therefore,  the  dimensions  of  the  mesh  opening  will  be  chosen  to  be  at  most  .1  A.  by 
.  1  A  . 

The  simulator  is  an  open  waveguide  structure,  and  consists  of  flat  conducting 
plates  that  are  finite  in  width.  As  observed  in  the  PEC  plate  problem,  a  large  amount 
of  current  will  be  distributed  near  the  edge  of  the  plate.  Therefore,  in  our 
configuration,  it  is  extremely  important  to  include  an  edge  wire  to  approximate  this 
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singular  current.  It  should  also  be  recognized  that  there  are  currents  on  both  sides  of 
the  conducting  surface  in  the  physical  problem.  However,  with  the  use  of  a  thin-wire 
mesh  approximation  the  induced  current  lies  on  the  axis  of  each  wire  segment. 
Furthermore,  the  current  on  each  segment  is  limited  to  one  dimension.  This  problem 
also  consists  of  a  symmetrical  structure  about  the  center  longitudinal  axis.  By 
imposing  symmetry  on  the  problem,  the  computation  time  of  the  matrix  inversion  can 
be  reduced  by  a  factor  of  N1! 4.  However,  in  doing  so,  an  odd  mode  solution  is 
enforced  in  the  transverse  direction,  and  Ex  is  zero  in  the  x-0  plane. 

The  source  and  load  configurations  will  be  the  same  as  those  developed  for  the 
rhombic  EMP  simulator  in  Chapter  2.  The  source  consists  of  a  constant  E-field 
distributed  over  wire  vertical  wire  segments  placed  symmetrically  about  the  y-axis. 
It  must  be  realized  that  the  desired  source  of  the  parallel-plate  simulator  is  an  E-field 
distributed  evenly  over  the  entire  source  gap  region. 

When  the  simulator  is  terminated  by  its  characteristic  impedance,  a  large  amount 
of  the  energy  available  is  dissipated  by  the  load  and  very  little  reflection  is 
encountered.  Since  the  load  is  to  be  matched  to  the  dominant  TEM  mode,  which 
contains  only  a  vertical  E-field  component,  the  loading  wire  segments  have  been 
placed  vertically  across  the  load  gap  region  in  a  similar  configuration  used  to 
terminate  the  rhombic  EMP  simulator.  Although  at  higher  frequencies,  where  TE  and 
TM  modes  will  propagate,  the  line  is  no  longer  matched,  and  less  energy  will  be 
dissipated  by  the  load.  As  a  result,  more  energy  has  the  potential  of  being  reflected 
back  into  the  simulator.  A  common  technique  used  to  minimize  reflections  is  to  use  a 
sloped  distributed  load.  In  doing  so,  not  only  the  vertical  E-field  is  matched,  but  the 
longitudinal  fields  of  the  TM  modes  will  be  matched  as  well.  Furthermore,  the  load  is 
generally  slightly  inductive  to  account  for  the  change  in  characteristic  impedance  at 
higher  frequencies.  Another  advantage  of  sloping  the  load  is  that  the  required  current 
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in  the  termination  will  reach  its  late  time  value  much  faster  during  pulse  excitation. 
This  can  be  explained  in  a  simplified  manner  by  the  fact  that  the  currents  entering  the 
load  from  the  conducting  plates  will  reach  a  position  on  the  terminating  load  much 
faster  after  the  arrival  of  the  incident  wave.  In  doing  so,  the  load  becomes  better 
matched  in  its  quasi-static  characteristic  impedance.  Furthermore,  since  the  mcidem 
wave  arrives  at  different  times  along  the  surface  of  the  load,  the  reflections 
encountered  are  dispersed  in  time,  and  the  incident  wave  is  reflected  back  at  an 
oblique  angle,  and  not  directly  back  into  the  simulator.  For  a  detailed  study  of  sloped 
terminations  of  parallel-plate  simulators,  the  reader  is  referred  to  [35],  [36],  and  [37], 
Our  analysis  with  NEC  will  not  be  able  to  consider  a  sloped  termination  due  to  the 
limitations  of  the  load  gap  model. 


3.2  Numerical  Results 

The  first  case  to  be  analyzed  is  that  of  an  asymmetrical  parallel-plate  simulator, 
with  conical  tapered  end  sections.  The  dimensions  will  be  modeled  after  the 
ACHATES  simulator,  which  was  designed  by  the  Air  Force  Weapons  Laboratory  [34], 
The  ACHATES  simulator  is  chosen  for  its  unique  characteristics,  which  will 
subsequently  enable  the  results  to  be  interpreted  easily. 

The  ACHATES  simulator  has  a  height  to  width  ratio  of 

1=1.2335.  (3.1) 

a 

This  aspect  ratio  is  chosen  such  that  the  simulator  has  a  quasi-static  characteristic 
impedance  of  100  Q.  Furthermore,  it  can  be  shown  that  one- volt  of  input  voltage 
across  the  source  gap  of  ACHATES  will  produce  a  field  of  one-volt  per  meter  at  the 
center  of  the  ground  plane  (at  x-O).  In  Chapter  1  the  bandwidth  of  a  parallel-plate 
simulator,  which  is  determined  by  the  dispersion  distance  associated  with  the  spherical 
wavefront  that  is  launched  into  the  parallel-plate  region  from  the  conical  tapered  end 


vs 


section,  wa  discussed.  It  was  found  that  the  more  gradual  the  taper,  the  higher  the 
bandwidth  of  the  simulator.  The  design  of  ACHATES  initially  considered  is  one  that 
has  very  long  tapered  end  sections,  resulting  in  a  bandwidth  of  450  MHz.  The 
dimensions  of  the  simulator  are  illustrated  in  Figure  3.1.  and  the  wire  mesh 
approximation  generated  for  NEC  is  illustrated  in  Figure  3.2. 

Three  different  operating  frequencies  are  chosen  for  this  example:  1 )  75  MHz.  2) 
100  MHz.  and  3)  125  MHz  excitations.  For  each  case  the  characteristic  impedance  of 
the  simulator  is  computed  by  determining  the  load  resistance  that  minimizes  the 

current  standing  wave  ratio  and  also  minimizes  the  radiation  efficiency.  Then  the 

r 

distribution1 0/  E-fields  in  the  center  parallel-plate  region,  which  is  often  referred  to  as 
the  "working  volume."  is  analyzed. 

Table  3.1  illustrates  the  characteristics  of  the  simulator  for  various  terminating 
loads.  The  computed  characteristic  impedance  is  120  Q  for  a  75  MHz  excitation. 
120  Q  for  a  100  MHz  excitation,  and  95  Q  for  a  125  MHz  excitation.  In  comparison  to 
the  rhombic  EMP  simulator,  the  characteristic  impedance  seems  to  be  much  more 
stable  over  a  broad  frequency  range.  Furthermore,  if  the  simulator  is  matched  into  a 
100  Q  load,  at  most,  an  SWR  of  1.3  in  the  current  will  be  experienced  over  the 
f  requency  range  considered. 

The  E-field  distribution  in  the  working  volume  due  to  a  one- volt  continuous 
wave  excitation  across  the  source  gap  is  illustrated  in  Figures  3.3  through  3.11.  The 
vertical  E-field,  Ex ,  along  the  y-direction,  is  illustrated  in  Figures  3.3  to  3.5  for  the 
three  cases  considered.  The  fields  are  quite  uniform  and  have  a  linear  progressive 
phase  shift  throughout  the  working  volume.  It  should  also  be  noted  that  the  phase 
has  little  change  along  the  vertical  direction  since  the  propagating  waves  have  planar 
equiphase  fronts.  Recall  that  the  design  of  ACHATES  is  such  that  one-volt  of  input 
voltage  produced  one-volt  per  meter  at  the  ground  plane  (at  x-0),  in  the  quasi-static 


case.  This  is  true  for  ail  three  cases.  The  deviation  is  due  to  the  existence  of  non-TEM 
modes,  which  become  more  significant  as  the  frequency  is  increased. 

Figures  3.6  to  3.8  illustrate  the  vertical  E-field  as  a  function  of  the  vertical 
direction.  The  fields  are  quite  uniform  throughout  the  height  of  the  simulator. 
However,  it  is  apparent  that  there  is  error  in  the  computed  field  close  to  the  thin-wire 
mesh.  Figures  3.9  to  3.11  illustrate  the  longitudinal  component.  Ey ,  in  the  x-0  plane 
and  is  representative  of  the  propagating  TM  field.  The  longitudinal  field  is  on  the 
order  of  a  magnitude  smaller  than  the  vertical  E-field.  However,  it  becomes  more 
significant  in  the  higher-frequency  regions.  This  can  be  explained  as  follows.  A 
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"quasi-spherical"  wavefront  propagates  from  the  wave  lauching  conical  tapered  end 
section  into  the  parallel-plate  region.  Along  the  ground  plane,  the  incident  E-field  is 
vertical  and  there  is  no  longitudinal  component.  Moving  above  the  ground  plane,  in 
the  x-0  plane,  the  field  consists  of  both  TM  and  TEM  modes.  Then,  getting  closer  to 
the  top  plate,  the  TM  component  becomes  more  significant  since  Ey  is  larger,  due  to 
the  characteristics  of  the  incident  wavefront.  Therefore,  at  higher  frequencies  the 
trailing  edge  of  the  spherical  wavefront  will  become  comparable  to  a  wavelength  and 
significant  TM  modes  are  generated.  This  happens  at  the  upper-frequency  cutoff  of 
the  simulator,  which  is  450  MHz  for  this  example.  Furthermore,  it  is  noted  that  £y 
tends  to  be  a  minimum  at  y43  and  increases  again  as  y  progresses  due  to  the  mismatch 
with  the  terminating  conical  end  section. 

In  order  to  develop  a  better  understanding  for  the  mechanism  contributing  to  the 
large  aberration  in  the  planar  wavefront  of  the  parallel-plate  region,  the  field 
distribution  within  the  simulator  is  analyzed  near  its  cutoff.  By  shortening  the 
conical  tapered  end  sections  from  4.5  m  to  1.95  m.  the  bandwidth  of  the  simulator  is 
reduced  from  450  MHz  to  200  MHz  [34).  In  doing  so,  the  simulator  has  a  computed 
characteristic  impedance  of  UOO  for  a  100  MHz  excitation  and  a  corresponding 


current  SWR  of  1.125.  Furthermore,  at  175  MHz.  which  is  very  near  cutoff,  the 
simulator  is  matched  when  loaded  by  95  Cl.  and  has  a  corresponding  current  SWR  of 
1 .08.  It  should  also  be  noted  that  the  radiation  efficiency  for  the  1 75  MHz  case  is  80. 19 
%.  Therefore,  verv  little  energy  is  being  dissipated  by  the  load  resistance,  and  most  of 
the  energy  is  either  being  reflected  or  radiated.  Figures  3.12  to  3.17  illustrate  the  E- 
field  distribution  in  the  working  volume  due  to  a  one-volt  excitation  across  the  source 
gap. 

The  vertical  E-field  as  a  function  of  y  is  illustrated  in  Figures  3.12  and  3.13.  At 
100  MHz.  the  field  is  uniform  and  has  a  linear  progressive  phase  shift.  However,  close 
to  the  simulator  cutoff  frequency,  there  is  a  large  contribution  from  higher-order  TM 
and  continuous  modes.  The  vertical  field  deviates  significantly  from  its  uniformity 
and.  in  fact,  has  a  1/R  contribution,  where  R  is  the  distance  from  the  source. 
Furthermore,  the  phase  shift  is  no  longer  linear.  It  should  be  noted,  although  it  is  not 
illustrated,  that  off  of  the  x-0  plane,  there  is  a  large  contribution  to  the  field  from 
higher-order  TE  modes  as  well.  The  distribution  of  the  vertical  E-field  in  the  vertical 
direction  is  illustrated  in  Figures  3.14  and  3.15.  At  100  MHz.  the  field  is  relatively 
uniform,  and  the  phase  is  approximately  constant  over  the  height  of  the  parallel  plate. 
However  at  175  MHz.  the  field  is  no  longer  uniform.  Furthermore,  the  phase 
indicates  that  the  modes  contributing  to  the  spherical  wavefront  are  propagating  with 
little  attenuation  into  the  working  volume.  At  y  equals  -0.25  m.  there  is  a  45  degree 
phase  lag  from  z  equals  0  to  z  equals  0.94  m.  Doing  a  rough  computation,  the 
maximum  deviation  of  the  spherical  wavefront  at  the  planar  interface  is  2.165  -  1.95. 
or  0.215  m  (where  2.165  is  the  length  the  wavefront  has  to  travel  along  the  top  plate 

of  'he  conical  taper).  At  175  MHz.  this  corresponds  to  a  x  .215  phase  shift  of 

A. 


45. 15  degrees.  The  wavef ront  becomes  planar  as  the  wave  propagates  with  y,  although 
it  loses  a  considerable  amount  of  energy. 
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Figures  3.16  and  3.17  illustrate  the  longitudinal  E-field  along  the  vertical 
direction.  Initially  comparing  Figure  3.16  to  Figure  3.10.  the  longitudinal  component 
of  the  E-field  is  much  larger  towards  the  waveguide  interface  at  y— .50  m.  However, 
at  175  MH/.  there  is  a  significant  increase  in  the  amplitude  of  Ey  throughout  the 
parallel-plate  region.  An  interesting  observation  in  Figure  3.17  is  that  the  location  of 
the  peak  amplitude  of  Ey  along  the  vertical  occurs  at  shorter  heights  as  the  wavefront 
travels  f rom  one  end  of  the  parallel-plate  region  to  the  other.  This  can  be  explained 
from  ray  optics.  The  incident  field  is  a  diverging  wavefront,  bounded  by  the  conical 
tapered  end  section.  Near  the  ground  plane,  the  waves  propagate  directly  into  the 
parallel-plate  region.  However,  for  larger  angles  of  incidence,  the  waves  are  reflected 
by  the  surface  of  the  top  parallel  plate.  The  maximum  angle  is  determined  by  the 
angle  of  the  tapered  end  section.  Once  reflected,  this  ray  determines  the  reflected 
shadow  boundary,  which  is  propagating  in  the  negative  z  direction  as  it  progresses  in 
y.  Therefore,  due  to  the  constructive  interference  of  he  longitudinal  component  of 
the  E-f leld.  the  peak  field  progressively  occurs  at  shorter  distances  above  ground. 

3.2.1  Computational  considerations 

All  computations  were  done  on  an  AT&T  PC6300  desktop  microcomputer  or  a 
compatible  DOS  machine.  The  NEC  code  was  rewritten  and  compiled  in  double 
precision  Microsoft  Fortran.  The  program  itself  requires  450  kbytes  of  RAM  to  load, 
and  approximately  8  Mbytes  free  on  a  hard  disk  for  out  of  core  solutions  of  matrices 
with  600  xJbOO  elements.  Table  3.2  lists  the  computation  times  required  to  fill  and 
factor  the  moment  matrix  for  the  ACHATES  simulator  with  a  450  MHz  bandwidth. 
The  computation  times  were  recorded  from  runs  on  the  AT&T  and  are  listed  in  Table 
1.2.  It  should  be  noted  that  the  above  computations  include  one  plane  of  symmetry 

V2 

about  the  x-axis,  which  reduces  the  computation  time  by  about  — — .  As  the 
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frequency  becomes  larger,  the  number  of  unknowns  becomes  significantly  larger,  and 
the  problem  becomes  extremely  time-consuming  to  solve.  Furthermore,  note  that  it 
will  require  3.5  Mbytes  to  store  a  matrix  of  466  by  466  in  double  precision.  In  the 

i 

high-frequency  region,  such  that  the  wavelength  is  shorter  than  the  parallel-plate 
|  height,  the  problem  becomes  too  big  to  be  practically  handled  by  NEC.  For  instance,  if 

we  increase  the  frequency  to  150  MHz,  this  will  require  656  unknowns,  nearly  7 
Mbytes  to  store,  approximately  22,982  seconds  to  fill,  and  approximately  33.b40 
seconds  to  factor. 

3.3  Interpolation  of  the  Wire  Mesh  Currents 

While  studying  the  PEC  plate  problem  in  Chapter  2.  it  was  found  that 
discontinuities  appear  in  the  wire  segment  currents  at  multiwire  junctions.  This  is  a 
result  of  enforcing  Kirchoff's  Law  as  a  condition  in  solving  for  the  unknown 
weighting  coefficients.  However,  a  discontinuity  in  the  current  will  lead  to 
nonphysical  singular  behavior  in  the  near  fields,  which  will  be  referred  to  as  “singular 
distortion."  Due  to  the  wire  mesh  approximation  of  the  conducting  surface  of  the  EMP 
simulator,  there  are  discontinuities  in  the  wire  current  as  well.  Figure  3.18  illustrates 
the  real  and  imaginary  components  of  the  current  along  the  center  wire  of  the 
simulator.  Each  step  discontinuity  is  located  at  a  multiple  wire  junction  in  the  wire 
mesh.  Figure  3.19  illustrates  the  resultant  scattered  E-field  computed  in  the 
tangential  direction  on  the  center  wire  surface  along  the  conical  wave  launching  end 
section.  Note  that  the  singularity  in  the  E-field  is  the  largest  at  the  source  gap.  This  is 
due  to  boundary  condition  error  suffered  by  the  source,  as  well  as  the  oblique  angle 
the  source  segment  makes  with  the  conical  tapered  end  section.  It  is  obvious  that  these 
singularities  are  nonphysical;  however,  they  are  a  trade-off  that  one  must  take  when 
using  NEC  to  solve  wire  mesh  scattering  problems.  In  Section  3.2,  it  was  obvious  that 
the  results  were  plagued  by  error  as  the  position  of  the  near-field  computation  moved 


closer  to  the  conducting  surface.  These  singularities  can  be  removed  by  smoothing  the 


current.  Therefore,  once  the  currents  are  smoothed,  it  can  be  determined  how  close  to 


the  wire  mesh  the  field  can  be  computed  before  suffering  significant  error  from  the 


leld  singularities. 


3.3.1  One-dimensional  smoothing  of  the  wire  currents 


Consider  a  one-dimensional  cunent  described  by  some  function  f(x)  whose 
support  lies  on  the  interval  (0.1. ).  The  function  f( x )  can  be  expanded  into  a  Fourier 


sme  series  as 


/  (. r  )  =  £  bk  sin. 


(3.3.1 


However,  the  expansion  of  f(x)  will  have  a  very  poor  convergence  since  the  boundary 


conditions  at  x-Q  or  x-L  will  not  be  fulfilled  unless  f(x)  vanishes  smoothly  at  the 


endpoints.  Next,  a  function  g(x)  is  introduced  such  that 


g  (0)  =  0 


(3.3.2) 


g(L  )  =  0 


g(-x  )  =  -g (  t  )  . 


Therefore.  g(x)  is  chosen  to  be 


g  (x  )  =  f  (x  )  -  (a  +bx  ) 


(3.3.3) 


where 


a  =  /  (0) 


b  =  f  (L  )-f  (0) 

As  a  result,  a  f unction  that  has  no  discontinuities,  and  is  periodic  with  a  period  of  2L. 


is  obtained.  A  discontinuity  does  appear  in  the  second  derivative  at  the  end  points; 


however,  any  singularities  in  the  first  two  derivatives  have  been  eliminated.  g(x)  can 


then  be  expessed  by  a  Fourier  sine  series  expansion  with  excellent  convergence  for  a 


well-behaved  fix).  Subsequently.  f( x )  can  be  obtained  f rom  g(x)  as 

f  [x  )  =  ?(  r  )  +  a  +  bx  .  (3.3.4 

The  next  step  is  to  discretize  fix)  at  equidistant  points,  with  a  separation  h 
Therefore,  choosing  N  discrete  points  on  the  interval  ( 0. L )  leads  to 


yt.  =  f  ikh  )  Ik  =0,1 . N  ) 


(3.3.5: 


yk  '  =  g  (kb  )  {k  =0.1. 
The  Fourier  series  is  expressed  as 


y*  '  =  L  hosin  I  "'V  '-  •  (3.3.7) 

i> = i  (  L 

The  constant  coefficients,  b  a.  are  solved  by  taking  the  inner  product  of  Equation 


a  nkh 
~L~ 


(3.3.7)  with  smfL-^L.  However,  since  g(kh)  has  a  discrete  finite  support  on  the 

interval  A:  €  [O./V  ].  and  the  constant  coefficients  can  be  derived  by  the  finite 
summation 


?  v 

bm  =  TT  L  V*  Sin 


(3.3.8) 


/v*=i  (  'v  j 

However,  the  spacing,  h.  of  the  discrete  points  and  the  number  of  coefficients  needed 
for  proper  convergence  must  be  established.  The  Fourier  series  expansion  is  a 
superposition  of  Fourier  harmonics.  There  exists  a  cutoff  frequency.  u0,  where  the 
function  g(x)  has  no  more  existing  harmonics  [38].  Therefore,  in  the  absence  of  noise. 

the  Fourier  coefficients.  bvb2 . by.  will  become  increasingly  smaller  above  some 

value  Hewever.  in  the  presence  of  noise,  the  higher-order  coefficients  are  no  longer 
negligible.  If  the  noise  were  random,  the  frequency  spectrum  would  be  a  constant 
over  all  frequencies.  Therefore,  the  higher-order  coefficients  are  relatively  constant  in 
magnitude  and  have  a  random  phase  which  introduces  undesirable  harmonic  distortion 
in  the  interpolated  function.  Therefore,  the  series  has  to  be  truncated  at  some  finite 
value  of  k  =^.  This  number  is  best  determined  by  analyzing  the  data  and  subjectively 
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choosing  the  cutoff  f  requency.  v0,  which  is  expressed  as  [38] 


n  Ji  _ ~ 

—  X  “2jriyc 


(3.3.0) 


m  , 

-  = ^ 


L  =  Mh  . 


The  real  and  imaginary  currents  to  be  interpolated  are  periodic  in  nature,  and  the 


noise  is  introduced  by  the  discontinuities  at  each  multiple  wire  junction.  Since  the 


function  is  relatively  well  behaved,  the  coefficients  diminish  rather  quickly. 


However,  due  to  the  noise,  the  harmonic  distortion  introduces  a  ringing.  By 


experimenting,  it  was  found  that  the  cutoff  frequency  is  best  chosen  to  be 


=  IT 


(3.3.10) 


This  leads  to  £  =  12  coefficients  per  wavelength. 


Figure  3.20  illustrates  the  interpolated  current,  in  the  longitudinal  direction,  for 


the  case  of  a  100  MHz  excitation  of  the  450  MHz  bandwidth  ACHATES  simulator. 


The  real  and  imaginary  currents  are  interpolated  separately  in  order  to  preserve  the 


phase.  It  is  easily  seen  that  smoothing  by  Fourier  sine-series  interpolation  provides 


excellent  results  for  this  case. 


The  transverse  currents,  which  are  on  an  order  of  magnitude  smaller  than  the 


longitudinal  currents,  are  subject  to  a  large  amount  of  relative  error.  This  was 


ibserved  in  the  PEC  patch  problem  discussed  in  Chapter  2.  Since  Kirchoff's  law  is 


being  satisfied  at  the  wire  junction,  small  discontinuities  m  the  longitudinal  current 


lead  to  relatively  large  discontinuities  in  the  transverse  current.  There  are  also  a  few 


other  limitations  that  are  encountered  due  to  the  wire  mesh  approximation.  The  first 


is  that  the  edge  wire  is  not  solely  in  the  y-direction  along  the  tapered  end  sections,  and 


it  has  a  vector  component  in  the  transverse  direction.  However,  the  vector 


.ontnbution  of  current  is  not  accounted  for  at  junctions  with  the  transverse  wires. 
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Furthermore,  the  current  does  not  vanish  in  the  x-0  piane,  due  to  the  discontinuity  of 
the  junction  of  the  transverse  wires  with  the  center  longitudinal  wire.  When 
interpolating  the  transverse  current  these  factors  will  have  to  be  taken  into  account. 


Recall  that  the  impedance  matrix  ot  the  moment  method  solution  was  generated 
by  point  matching.  Therefore,  the  boundary  condition  on  the  tangential  E-field  is 
being  satisfied  at  the  center  ot  each  wire  segment.  Furthermore,  the  large 
discontinuity  being  suffered  at  either  endpoint  is  due  to  the  perpendicular  longitudinal 
currents.  As  a  result,  only  the  segment  center  currents  are  interpolated  while 
enforcing  the  end-point  boundary  conditions  and  imposing  odd-symmetry  across  the 
x-0  plane.  Finally,  the  real  and  imaginary  currents  are  interpolated  separately  by  a 
Fourier  sine-series  expansion  as  was  previously  described.  Figures  3.21  and  3.22 
illustrate  the  interpolated  transverse  current  magnitude. 

Once  the  currents  are  interpolated,  the  near  E-field  can  be  easily  computed  with 
the  use  of  the  integral  differential  equation  derived  from  Maxwell's  equations. 
However,  since  the  currents  are  one-dimensional,  the  current  and  the  Green's  function 
are  parameterized,  leaving  only  one  variable  of  integration.  The  resultant  equation  is 
given  as 

T 

E{r)=-±^tr  f  ( VV-f"  +k2)I  (t  )G(r  .t  ')di  (3.3.11) 

4t rkz  Jo 

where 


C(c  .;  )  = 


R  =  ((.t  —  x  )2+( y  —  y  )2+(r 


v  =  x0l  +al 

y '  =  yot  +to 

2  '  =  -o. 

(  xot ,  v,,  .  z0l  )  is  the  position  of  the  wire  starting  point  in  cartesian  coordinates,  and 


p— p'  =  x  { x  —x  ' )  +  y  ( y  —y  )  +  z  i  z  —z  ’ )  . 

Figure  3.22  illustrates  an  example  of  the  E-field  computed  from  the  interpolated 
currents.  Initially,  comparing  Figure  3.22  to  Figure  3.4,  it  is  found  that  singular 

distortion  in  the  vertical  E-field  is  negligible  at  heights  up  to  where  h  is  the  height 

of  the  parallel-plate  above  the  ground  plane.  At  a  height  of  .  it  appears  that  the 

4 

singular  distortion  in  Figure  3.4  is  still  tolerable.  Furthermore,  it  is  interesting  to  note 
that  there  is  very  little  distortion  of  the  phase  due  to  the  presence  of  the 
discontinuities  of  the  current. 


It  must  be  realized  that  interpolating  the  currents  and  enforcing  the  edge 
conditions  do  not  improve  the  boundary  condition  error.  In  fact,  it  only  eliminates  the 
singulantiesmn  the  near  electric  field.  However,  this  analysis  does  give  insight  into 
how  much  the  singularities  contribute  to  the  fields  when  computed  directly  by  NEC. 
Therefore,  combining  these  results  with  those  in  Section  3.2,  it  is  found  that  a  useful 
and  valid  approximation  of  the  field  distribution  in  the  working  volume  has  been 
provided  by  NEC.  However,  singular  distortion  occurs  when  the  fields  are  computed 


near  the  wire  mesh.  In  the  previous  example,  singular  distortion  has  a  significant 


contribution  when  the  fields  are  computed  within  _  ot  the  top  plate  (whet\  he 
height  of  the  simulator  is  approximately  -i-  at  100  MHz).  Therefore,  when  applying  a 

thm-wire  mesh  approximation  to  a  stir! ace  scattering  problem,  one  should  be  aware  ol 
the  singular  distortion  before  interpreting  the  results. 


j.4  Simulator-Obstacle  Interaction 


The  purpose  of  the  EMP  simulator  is  to  reproduce  the  electromagnetic 
environment  of  the  EMP  radiated  by  a  high  altitude  nuclear  burst  and  to  subject 
various  test  objects  to  this  environment.  Thus  far.  the  characteristics  of  the  parallel- 
plate  bounded  wave  EMP  simulator  have  been  discussed.  The  fields  within  the 
working  volume  of  the  simualtor  approximate  those  of  the  free-space  EMP 
environment.  However,  once  a  test  object  is  placed  within  the  working  volume  of  the 
simulator  the  characteristics  of  the  simulator  may  change.  From  a  quasi-static 
analysis  it  can  be  shown  that  the  object  capacitively  couples  to  the  upper  conducting 
plate  and  the  ground  plane.  Since  the  object  is  of  finite  length,  it  acts  as  a  perturbation 
of  the  transmission  line.  Dynamically,  the  obstacle  is  a  scattered  and  the  fields  are  a 
superposition  of  incident  fields  and  scattered  fields,  although  multiple  scattering 
occurs  between  the  simulator  and  the  obstacle.  In  both  analyses  it  is  this  interaction 
between  the  obstacle  and  the  simulator  that  causes  a  distortion  ol  the  simulator 
characteristics.  It  is  the  purpose  of  this  section  to  investigate  the  changes  in  the 
characteristics  of  the  EMP  simulator  when  placing  a  large  object  in  the  working 
volume. 

The  obstacle,  or  test  object,  is  a  cube  with  a  side  length  of  h/2  (where  h  is  the 
height  of  the  parallel  plate  above  ground).  The  cube  will  be  placed  symmetrically  in 
the  parallel-plate  region  between  the  top  plate  and  ground,  and  in  the  x-direction. 
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Furthermore,  the  front  face  of  the  cube  will  be  placed  in  the  transition  plane  ol  the 
conical  tapered  end  section  and  the  parallel-plate  region.  This  is  illustrated  in  Figure 
3.23.  As  an  example,  the  450  MHz  bandwidth  ACHATES  EMP  simulator,  excited  by  a 
100  MU/  ".mi  amplitude  excitation  is  used. 

For  small  simulator-obstacle  interaction,  negligible  changes  in  the  characteristics 
of  the  simulator  are  expected,  such  as  the  characteristic  impedance,  input  impedance, 
or  the  radiation  efficiency.  However,  in  the  presence  of  the  object,  a  perturbation  ol 
the  near  fields  is  still  expected  in  the  working  volume  since  the  fields  are  a 
superposition  ol  the  incident  field  and  the  scattered  fields.  Referring  to  Table  3.1.  the 
matched  load  condition  for  a  100  MHz  excitation  of  the  450  MHz  bandwidth  simulator 
with  an  empty  working  volume  is  RL  =  120  ft.  Furthermore,  the  simulator  has  an 
input  impedance  of  97.04+ j  14.68  ft  and  a  radiation  efficiency  of  20.77  %.  Once  the 
cube  is  placed  in  the  working  volume,  as  in  Figure  3.23.  the  input  impedance  reduces 
to  76.64+jJ  7.03  ft.  By  lowering  the  impedance,  the  available  input  power  is  increased, 
however,  it  is  interesting  that  the  radiation  efficiency  reduces  to  Id.  14  %.  This 
indicates  that  a  larger  amount  of  power  is  being  dissipated  by  the  load  impedance. 
Figure  3.24  illustrates  the  longitudinal  currents  on  the  top  plate  of  the  simulator  for 
the  case  of  an  empty  working  volume,  and  when  the  obstacle  is  in  the  working 
volume.  Comparing  Figure  3.24(a)  with  Figure  3.24(b)  it  is  observed  that  a  slight 
change  in  the  longitudinal  current  is  encountered  with  the  presence  of  the  object.  This 
slight  change  in  currents  introduces  a  small  change  in  the  held  distribution 
characteristics  of  the  simulator,  which  can  essentially  be  referred  to  as  the  incident 
field.  Therefore,  the  field  within  the  simulator  is  a  superposition  of  this  incident  field 
and  the  scattered  field  produced  by  the  currents  induced  on  the  test  object.  Figures 
T25  to  T27  illustrate  the  comparison  of  the  field  distributions  lor  the  cases  of  an 
empty  working  volume,  or  free-space  environment,  and  in  the  presence  of  an  obstacle, 
or  a  perturbed  f teld  environment. 


Following  the  above  analysis,  there  is  additional  work  that  could  be  done  to  help 
characterize  simulator-obstacle  interaction.  In  particular,  studying  the  changes  in  the 
simulator  characteristics  over  a  broad  frequency  range  is  necessary  to  interpret  the 
uroad  uand  excitation  ot  the  test  object.  Furthermore,  computation  ot  the  fields  due 
*o  the  perturbed  simulator  currents  and  in  the  absence  of  ’he  obstacle  currents  could 
provide  some  useful  results  in  determining  the  perturbation  of  the  incident  field  in  the 
presence  of  the  obstacle.  Finally,  an  analysis  of  the  plane-wave  scattering  of  the 
obstacle,  in  the  absence  of  the  top  conducting  plate,  could  provide  useful  results  as 


TABLE  3.1(c) 

NEC  MODEL  OF  THE  ACHATES  EMP  SIMULATOR:  125  MHz 
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1.837 
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56.81 

3.504 

1000.0 

43.67 

55.02 

71.86 

5.115 

5000.0 

35.84 

57.66 

92.16 

8.380 

TABLE  3.2  COMPUTATIONAL  CHARACTERISTICS  OF  NEC 


( run  on  an  AT&T  PC6300  microcomputer) 
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Interpolated  vs.  NIC  transverse  Current 
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l-'igiirt'  <  21.  Comparison  of  the  currents  on  the  transverse  wires  that  are  smoothed  by  a  I  ourier  sine 
series  mterpolation.U'urrenls  are  normalized  by  a  wavelength) 


Imaginary  part,  at  y--0.5  m,  with  the  edge  condition  imposed  by  the  edge  wire. 
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in  the  presence  of  the  obstacle. 
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CHAPTER  4 


CONCLUDING  REMARKS 

I  he  p  '  "a  t.'-m''  w  ds  »o  develop  the  characteristics  of  a  parallel-plate 

''•'•incit'd  r  1 1  a i o r  by  scattering  theory,  using  a  thin-wire  mesh 

appmximah  re  '.inducting  surfaces  and  the  Numerical  Electromagnetics  Code 

(NEC)  [22].  1’  'v as  found  that  NEC.  which  is  a  user-oriented  computer  code  that 
offers  a  moment  me’ hod  solution  of  thin-wire  scattering  problems,  can  be  an 
extremely  useful  tool  for  analyzing  this  problem.  In  order  to  develop  the  wire  mesh 
model  and  an  understanding  of  the  limitations  associated  with  it.  Chapter  2  focused  on 
developing  corn  plicated  thm-wire  scattering  problems  that  in  turn  could  be  related  to 
the  parallel-pla’e  simulator.  The  emphasis  was  placed  on  determining  the  limitations 
of  the  wire  mesh  approximation  and  developing  a  proper  source  and  load  model  for  the 
parallel-plate  simulator. 

In  Chapter  2,  the  rhombic  EMP  simulator  was  analyzed.  The  rhombic  simulator, 
like  the  parallei-plate  simulator,  is  a  transmission-line  type  of  simulator.  The  source 
and  load  models  were  developed  using  a  thin-wire  approximation,  and  the  resultant 
haracteristics  of  the  simulator  were  compared  to  theoretical  and  experimental  data 
published  by  Shea  and  King  of  Harvard  University  [28,29].  The  general 
characteristics,  as  well  as  the  thin-wire  conducting  current  and  electric  field 
dis‘ n buttons  within  the  center  region,  compared  extremely  well  with  Shen  and  King  s. 
It  was  found  that  at  very  high  frequencies,  when  the  source  and  load  gaps  are  greater 
than  0. 1  a.  significant  error  in  the  value  of  the  matched  terminations  occurs  due  to  the 
wire  segment  approximation  of  the  source  and  load.  It  is  interesting  that  even  though 
'he  value  of  the  matched  termination  was  in  error,  the  currents  and  fields  still 
compared  extremely  well  with  the  results  of  Shen  and  King  under  matched  conditions. 


In  Chapter  2  the  thm-wire  mesh  approximation  of  a  square  PEC  plate  scatterer 
as  employed  m  order  to  evaluate  the  limitations  encountered  in  approximating  a 
'  ‘ e  *.  ;dt h .  lat  eond  icting  surtace.  I  he  plate  was  illuminated  by  a  normally 

,u  ’  If  mane  wave  and  ’bee  urrents  induced  on  the  plate  were  computed  bv  NEC 
n  ’  “e  w  ■  re  mesh  approximation.  I  hese  results  were  compared  to  those  derived  by  a 
'■  pec ;  r a i  Calerkin  ’echmque.  solved  by  conjugate  gradients.  I  he  wire  mesh 
approximation  provided  an  excellent  approximation  of  the  surface  currents.  However 
’  u.d  not  reproduce  the  edge  singularity  even  when  a  reasonable  number  of  unknowns 
were  used  to  represent  the  current. 

I  he  incident  E-lield  was  chosen  to  be  polarized  in  a  direction  parallel  to  an  edge 
>!  'he  plate,  which  is  referred  to  as  the  longitudinal  direction.  This  was  done  such 
'hat  the  currents  would  be  highly  oriented  in  one  direction,  and  extremely  small  in 
'he  transverse  direction.  The  transverse  curry's  vhich  are  at  least  an  order  of 
magnitude  smaller  than  the  longitudinal  currents,  suffered  a  large  amount  of  error. 
The  reason  for  the  error  is  due  to  the  limited  number  of  unknowns  as  well  as  large 
discontinuities  encountered  at  the  multiple  wire  junctions  in  the  wire  mesh. 
Furthermore,  the  discontinuities  are  a  result  of  the  enforcement  of  Kirchoff’s  current 
law  on  the  axial  currents  entering  each  multiple  wire  junction.  These  current 
discontinuities,  which  exist  in  both  the  transverse  and  longitudinal  directions,  will 
lead  to  nonphysical  singularities  in  the  near  electric  field.  Therefore,  one  must  be 
Lareful  as  to  how  close  the  E-field  can  be  computed  to  the  wire  mesh  surface  before 
suffering  from  singular  distortion. 

In  Chapter  3,  the  conducting  surface  of  the  parallel-plate  bounded  wave 
simulator  was  approximated  by  a  wire  mesh  which  was  developed  using  the  results 
obtained  from  the  examples  studied  in  Chapter  2.  The  well-defined  geometry  of  the 
ACHATES  EVIP  simulator  (34)  was  modeled.  T1  e  NEC  was  used  to  solve  tor  the 


currents  induced  on  the  wire  mesh  due  to  a  continuous  wave  excitation  across  the 
source  gap.  Using  these  results,  the  characteristics  of  the  simulator  were  established 
'"•er  a  oroad-t  requencv  region.  I  he  I  it* Id  distributions  in  ’he  parallel-plate  region 
under  mauled  or.d.tions.  were  computed  and  shown  to  be  predominately  FFM  tor 
.ow  t  req uencies  with  a  linear  progressive  phase  shift.  However,  at  higher  frequencies, 
w  hen  -he  operating  wavelength  was  near  that  of  the  cutoff  frequency  of  the 
simulator,  higher-order  modes  became  significant  due  to  the  mismatch  of  the  conical 
'apered  end  section  and  the  parallel-plate  waveguide.  As  a  result,  the  field 
distribution  greatly  deviates  from  the  desired  uniform  charactenst.es.  Finally,  the 

fr 

wire  mesh  model  was  used  to  analyze  the  changes  in  the  characteristics  of  the 
simulator  when  a  large  test  object  was  placed  in  the  working  volume.  It  was  found 
that  a  cube  of  side  length  equal  to  one-half  of  the  height  of  the  parallel  plate  above 
ground,  placed  symmetrically  in  the  working  volume,  has  a  small  effect  on  the 
characteristic  impedance  and  the  surface  currents  of  the  simulator.  This  small  change 
introduces  a  small  distortion  to  the  field  incident  upon  the  test  object.  The 
simulator-obstacle  interaction  is  important  to  understand  when  using  the  EMP 
simulator  as  a  testing  device,  and  is  worthy  of  further  investigation  using  similar 
'echniques  described  in  this  thesis. 

Overall,  the  Numerical  Electromagnetics  Code  provides  an  extremely  useful,  and 
relatively  simplistic,  analysis  of  the  parallel-plate  bounded-wave  simulator. 
However,  it  is  an  approximation  that  is  not  without  its  limitations.  Initially  the 
source  and  load  gaps  are  limited  to  a  one-dimensional  thin-wire  model  of  the  two- 
dimensional  gap.  For  the  rhombic  simulator,  this  had  an  effect  on  the  characteristic 
impedance  at  higher  frequencies.  Furthermore,  the  wire  mesh  approximation 
introduces  discontinuities  in  the  currents  which  lead  to  nonphysical  singularities  in 
the  near  electric  fields,  whereas  it  was  determined  that  fields  computed  at  distances 

beyond  JL.  from  the  surface  suffer  little  singular  distortion.  Finally,  the  number  ol 


unknowns  needed  to  approximate  the  conducting  currents  can  become  very  large,  since 
the  simulator  is  generally  much  larger  than  a  wavelength  in  the  mid-frequency  region. 
T  his  limits  the  analysis  ol  the  high-!  requency  region  tor  large  simulators. 
(  onsequerolv  despite  the  limitations  mentioned  above,  the  technique  ol  analysis 
•ne.n’ioned  herein  pros  ides  a  useful  'ool  in  the  analysis  of  the  parallel-plate  bounded  - 
wave  fcMP  simulator  in  the  low-  and  mid-frequency  regions. 
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